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Abstract 

In  this  research  the  electromagnetic  scattering  of  a  plane  wave  from  a  two-dimensional 
cavity  embedded  in  an  infinite,  perfectly  conducting  groimd  plane  is  investigated.  The 
plane  wave  is  assumed  to  be  xmder  transverse  electric  (TE)  polarization  with  respect  to  the 
s-axis.  The  cavity  may  be  empty  or  filled  with  an  arbitrary  homogeneous,  lossy  material. 
A  coupled  set  of  scalar  integral  equations  that  govern  the  electromagnetic  scattering  is 
implemented. 

An  approximate  solution  to  the  scalar  integral  equations  is  found  via  a  Method 
of  Moments  (MoM)  algorithm.  The  algorithm  is  implemented  in  a  computer  code,  and 
approximations  to  the  total  magnetic  field  on  the  cavity  smface  and  aperture  as  well  as  the 
normal  derivative  of  the  total  magnetic  field  on  the  cavity  aperture  are  obtained.  These 
fields  are  then  used  to  calculate  the  two-dimension^ll  monostatic  RCS  signatures  of  various 
test  cavities. 

The  numerical  results  from  the  algorithm  are  shown  to  agree  well  with  the  RCS 
signatures  calculated  by  other  well-known  methods  and  published  results.  In  addition  to 
being  accurate,  the  algorithm  is  very  computationally  efiicient.  The  process  results  in 
simply  solving  a  relatively  small,  well-conditioned  matrix  system  for  each  incident  angle 
to  produce  the  unknown  fields.  Analytical  and  numerical  error  analysis  shows  the  method 
to  have  approximately  0{h)  order  of  convergence,  where  h  is  an  indication  of  the  fineness 
of  the  mesh  used  to  discretize  the  cavity. 
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ANALYSIS  AND  NUMERICAL  SOLUTION  OF  AN  INTEGRAL  EQUATION 
METHOD  FOR  ELECTROMAGNETIC  SCATTERING  FROM  A  CAVITY  IN  A 

GROUND  PLANE 

/.  Introduction 

1.1  Motivation 

The  study  of  electromagnetic  plane  wave  scattering  and  the  measurement  of  the  radar 
cross  section  (RCS)  of  targets  has  been  of  much  interest  in  recent  years  to  the  military 
aircraft  industry.  Many  practical  applications  arise  which  require  knowledge  of  the  RCS 
signature  of  a  target.  The  RCS  of  a  target  is  essentially  a  measure  of  the  detectability  of 
the  target  by  a  radar  system.  More  rigorously,  in  three  dimensions,  the  RCS  or  echo  area 
is  defined  as  the  area  intercepting  the  amount  of  power  that,  when  scattered  isotropically, 
produces  at  the  receiver  a  density  that  is  equal  to  the  density  scattered  by  the  actual  target 
[6].  In  two  dimensions  the  RCS  parameter  is  referred  to  as  the  echo  width  or  alternatively 
as  the  radar  cross  section  per  unit  length  of  the  target. 

The  accurate  prediction  and  calculation  of  the  RCS  signature  of  targets  is  vitally 
important  to  today’s  military.  The  ability  to  control  the  RCS  signatures  of  friendly  air 
vehicles  as  well  as  the  ability  to  detect  and  identify  enemy  air  vehicles  are  both  significant 
determining  factors  in  the  outcome  of  an  air  campaign.  The  camoufiage  of  stealthy  aircraft 
is  one  example  of  the  need  to  minimize  a  friendly  target’s  RCS  signature.  Stealth  aircraft 
have  proven  to  be  a  highly  effective  means  of  utilizing  the  element  of  surprise  in  an  attack 
while  providing  increased  survivability  for  American  pilots,  as  evidenced  by  their  success 
in  the  Persian  Gulf  War. 

It  may  also  be  desirable  in  certain  instances  to  be  able  to  enhance  the  RCS  signature 
of  a  friendly  air  vehicle.  For  example,  the  use  of  unmanned,  remotely-piloted  air  vehicles 
to  gather  reconnaissance  information  and/or  serve  as  decoys  to  confuse  enemy  defenses 
may  be  a  very  effective  technique  in  certain  tactical  scenarios.  Such  a  vehicle  would 
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typically  be  much  smaller  than  a  fighter-size  aircraft  and  consequently  would  have  a  smaller 
RCS  signature.  The  ability  to  enhance  the  RCS  signature  of  a  smaller  unmanned  air 
vehicle,  thereby  reducing  an  enemy’s  capacity  for  distinguishing  between  it  and  a  fighter- 
size  aircraft,  may  be  advantageous. 

In  either  situation,  whether  reducing  the  RCS  or  enhancing  the  RCS  of  a  target, 
it  is  necessary  to  have  efiicient  methods  of  calculating  the  RCS  signature  of  a  scattering 
body.  Complex  objects  such  as  military  aircraft  have  many  different  design  factors  that 
can  contribute  to  their  RCS  signature,  but  of  particular  interest  is  the  prediction  of  the 
RCS  patterns  of  cavities  on  the  target’s  body.  Cavities  on  an  aircraft  can  include  jet 
engine  inlet  ducts  and  exhaust  nozzles,  cavity-backed  antennas,  as  well  as  cracks  and  gaps 
in  the  metallic  body  of  the  aircraft.  Cavities  are  of  such  importance  because  the  RCS 
of  a  cavity  can  potentially  dominate  the  total  RCS  signatme  of  an  aircraft,  and  because 
the  accurate  prediction  of  cavity  RCS  can  be  very  computationally  challenging.  Indeed, 
there  is  a  wealth  of  mathematical  and  engineering  literature  on  the  design  and  analysis  of 
electromagnetic  cavities.  (See,  for  example  [1],  [2],  [3],  [11],  [16],  and  [24].)  It  would  be 
far  too  expensive  and  time  consuming  to  rely  solely  on  real  world  measurements  in  the 
study  of  RCS  reduction  and  enhancement,  so  accurate  mathematical  models  and  efficient 
computational  algorithms  capable  of  predicting  RCS  patterns  are  much  desired. 

1.2  Organization 

The  goal  of  this  research  is  to  develop  an  efficient  and  accurate  method  of  calculating 
the  scattered  fields  of  an  incident  plane  wave  impinging  on  an  open  cavity  embedded  in 
a  grormd  plane.  These  scattered  fields  can  then  be  used  to  determine  the  RCS  of  the 
cavity.  To  accomplish  this  goal,  a  set  of  coupled  integral  equations  (developed  in  [26])  are 
implemented  and  solved.  In  Chapter  II,  the  cavity  geometry  and  the  problem  statement 
are  discussed  as  well  as  the  fundamental  concepts  from  electromagnetics  and  the  conditions 
satisfied  by  all  of  the  fields  in  the  problem.  Chapter  III  contains  the  development  of  the 
vector  integral  equations  that  govern  the  electromagnetic  scattering  of  a  plane  wave  from  a 
cavity.  The  vector  integral  equations  are  then  converted  to  a  set  of  scalar  integral  equations. 
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Chapter  IV  contains  a  brief  introduction  to  the  general  Method  of  Moments  (MoM) 
procedure  along  with  the  approximate  solution  of  the  sccdar  integral  equations  through 
the  development  of  a  MoM  algorithm.  Chapter  IV  also  contains  an  explanation  of  how 
certain  singular  integrals  arise  from  the  theory  and  how  they  can  be  treated  so  that  they 
evaluate  munerically.  An  estimate  of  the  order  of  convergence  of  the  MoM  algorithm  is 
also  presented  in  this  chapter.  The  niimerical  results  are  presented  in  Chapter  V,  including 
plots  of  the  fields  along  the  cavity  aperture  and  cavity  surface,  plots  of  the  monostatic 
RCS  signatures  of  several  test  cavities,  and  relative  error  plots  of  the  method.  Chapter  VI 
contains  a  summary  of  the  results  found  along  with  recommendations  for  future  research. 

1.3  Related  Work 

The  problem  of  predicting  the  electromagnetic  scattering  from  a  cavity  in  a  groxmd 
plane  has  been  fairly  well  studied,  and  there  are  many  different  techniques  and  approaches 
to  solving  the  problem.  As  expected,  these  techniques  have  their  advantages  as  well  as  their 
disadvantages.  Some  are  more  computationally  efiicient  then  others,  while  some  provide 
more  stable  and  accmate  solutions.  Presented  here  is  a  very  brief  overview  of  some  of  the 
more  popular  methods. 

The  generalized  network  formxilation  (GNF)  was  developed  by  Harrington  and  Mautz 
in  1976  to  calculate  the  scattering  from  an  open,  perfectly  conducting  cavity  in  free  space 
[14].  Their  method  partitions  the  cavity  problem  into  two  “simpler”  problems,  the  exterior 
cavity  problem  and  the  interior  cavity  problem.  Each  of  the  two  simpler  problems  is 
treated  individually  using  equivalence  relationships,  then  the  composite  problem  is  solved 
by  enforcing  the  continuity  of  the  tangential  fields  across  the  aperture  of  the  cavity.  The 
primary  drawback  to  the  GNF  method  is  that  the  fields  are  not  uniquely  determinable 
at  a  co;mtably  infinite  set  of  discrete,  resonant  frequencies  [21].  At  such  frequencies  a 
spurious,  nonsensical  solution  is  produced  by  the  method.  And  since  it  is  extremely  diflBcult 
to  determine  the  resonant  frequencies  of  an  arbitrarily  shaped  cavity,  this  is  a  rather 
significant  drawback.  Despite  this  hindrance,  however,  the  GNF  method  has  served  as  a 
springboard  for  many  other  methods.  In  fact,  most  integral  equation  approaches  employ 
some  form  of  a  GNF  implementation. 
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Other  populcu-  techniques  used  to  solve  the  scattering  problem  have  included  the 
finite  element  method  or  the  more  popular  hybrid  finite  element  methods.  These  hybrid 
methods  often  combine  traditional  finite  elements  with  boundary  integral  equations  to  solve 
the  scattering  problem  from  a  cavity  in  a  groimd  plane  [18].  They  have  the  advantage  of 
being  able  to  model  cavities  filled  with  arbitrary,  inhomogeneous  materials  since  the  entire 
cavity  interior  is  discretized  with  finite  elements.  Also,  the  spurious  resonance  problem 
discussed  earlier  can  be  avoided  by  introducing  penalty  functions  to  the  equations  [22]. 
However,  the  methods  are  typically  more  computationally  burdensome  due  to  the  required 
meshing  of  the  entire  cavity  interior  and  exterior.  But  this  is  eidmittedly  becoming  less 
and  less  of  a  problem  as  the  performance  of  computers  and  the  ease  of  the  storage  of 
information  are  dramatically  improving. 

Integral  equation  methods  represent  another  very  popular  family  of  techniques  for 
solving  the  cavity  scattering  problem.  The  electric  field  integral  equation  (EFIE)  and  mag¬ 
netic  field  integral  equation  (MFIE)  methods  are  easy  to  implement,  but  typically  only 
work  well  for  simple  geometries.  As  mentioned  earlier,  since  most  integral  equation  meth¬ 
ods  are  based  on  the  GNF,  they  suffer  from  the  problem  of  spurious  resonances.  However,  a 
new  integral  equation  method  was  introduced  by  Asvestas  and  Kleinman  [4]  in  1994.  They 
developed  a  set  of  coupled  vector  integral  equations  for  a  three-dimensional  imfilled  cavity- 
backed  aperture  in  a  perfectly  conducting  ground  plane,  which  were  purportedly  uniquely 
solvable  at  all  frequencies.  In  his  dissertation.  Wood  extended  Asvestas  and  Klenman’s 
coupled  vector  integral  equations  to  be  able  to  account  for  material  filled  cavities  [26]. 
The  method  he  developed  can  be  applied  to  electrically  large  and  electrically  small  cavi¬ 
ties,  and  it  will  not  suffer  from  the  problem  of  spurious  resonances.  Also,  the  technique 
requires  solving  only  for  the  imknown  fields  along  the  cavity  surface  and  apertxue,  so  it 
is  very  computationally  efficient.  Based  on  the  work  presented  in  his  dissertation.  Wood 
and  Wood  developed  a  related  set  of  scalar  integral  equations  for  the  scattering  from  a 
material  filled  two-dimensional  cavity  in  a  ground  plane  of  an  electromagnetic  wave  xmder 
transverse  magnetic  (TM)  polarization  [27].  This  current  research  effort  focuses  on  the 
complement  of  this  problem,  considering  an  incident  field  imder  transverse  electric  (TE) 
polarization. 
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II.  Theoretical  Background 

In  this  chapter,  the  general  geometry  of  the  embedded  cavity  is  presented  along  with 
the  formal  problem  statement.  Some  fundamental  electromagnetic  theory  is  then  given, 
followed  by  the  conditions  that  must  be  satisfied  by  all  of  the  fields  in  the  scattering 
problem.  Finally,  the  Green’s  functions  and  dyadic  functions  necessary  for  the  solution  of 
the  problem  are  defined. 

2.1  Cavity  Geometry 

Consider  the  geometry  of  the  two-dimensional  cavity  shown  in  Figure  1.  The  aper¬ 
ture,  or  cavity  opening,  is  denoted  as  <t,  and  the  smrface  of  the  cavity  is  called  S.  The 
interior  of  the  cavity  is  denoted  as  D  with  boundary  dD  =  crU5'  and  outward  pointing  imit 
normal  vector  n.  The  normal  vector  n  is  defined  almost  everywhere,  that  is  it  is  defined 
all  along  dD  except  for  a  finite  number  of  edges,  corners,  or  tips.  The  ground  plane  and 
the  surface  of  the  cavity  are  perfect  electric  conductors  (PEC),  and  the  upper  half  plane  is 
free  space  with  electric  permittivity  eo  and  magnetic  permeability  jUQ.  The  interior  of  the 
cavity  D  may  either  be  empty  or  filled  with  a  homogeneous  material  having  permittivity 
ei  and  permeability  ixy.  The  complement  of  a  in  the  ground  plane  is  called  cr*’. 
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The  symbol  r  will  be  used  as  a  position  vector,  that  is  r  =  xx  +  yy  +  2;z  in  three 
dimensions.  For  the  two-dimensional  problem,  the  projection  of  the  position  vector  into 
the  y^r-plane  is  used.  It  is  of  the  form  p  —  r-  (yy  -|-  zz)  —yy  +  zz.  The  image  of  an  object 
will  be  defined  as  the  reflection  of  the  object  across  the  groimd  plane  or  aperture.  For 
example,  the  images  of  S  and  D  in  Figure  1  are  Si  and  Dj,  respectively.  The  image  of  a 
point  given  hy  p  =  yy  +  zz  is  the  point  pi  =  p-ii  =  yy  —  zz,  where  Ij  =  yy  —  zz  is  the 
image  dyadic.  The  image  of  an  object  will  always  be  denoted  with  the  subscript  i  after  it. 
Similarly,  a  field  vector  function  is  defined  to  be  A(r)  =  Aj(r)x  -f  Ay(r)y  -I-  Az {r)z,  and 
its  image  is  defined  as  Ai(i^  =  Ax(r)x  +  Ay(f)y  —  Az{r)z. 

2.2  Problem  Statement 

In  the  problem  considered,  a  known  incident  plane  wave  whose  electric  and  magnetic 
fields  are  represented  by  impinges  on  the  cavity  giving  rise  to  the  unknown 

scattered  electric  and  magnetic  fields  (E®,  H®)  as  shown  in  Figure  2.  The  incident  plane 
wave  propagates  at  an  incident  angle  0  <  9  <  with  respect  to  the  positive  y-axis. 
Under  transverse  electric  (TE)  polarization,  the  magnetic  field  of  the  incident  plane  wave 
is  parallel  to  the  longitudinal  x-axis,  that  is  0).  The  electric  field  of  the 

incident  wave  has  both  a  y-  and  a  ^-component. 
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The  given  incident  fields  are  defined  to  be  those  fields  that  exist  in  unbounded  free 
space.  The  scattered  fields  are  those  that  are  caused  by  the  presence  of  the  scatterer,  in 
this  case  the  embedded  cavity.  They  are  defined  in  terms  of  the  total  field  and  the  reflected 
field  as  E®  =  E  -  E'”'^  -  E"®^  and  H®  =  H  -  H*"®  -  where  E  =  E(f)  is  the  total 
electric  field  intensity  in  the  region  in  volts  per  meter,  and  H  =  H(r)  is  the  total  magnetic 
field  intensity  in  the  region  in  Amperes  per  meter.  The  reflected  fields  E*^®^  and  H''®^  are 
defined  to  be  those  scattered  by  an  unbroken,  PEC  groimd  plane  located  at  =  0. 

2.3  Field  Conditions 

The  upper  half  plane  and  the  interior  of  the  cavity  are  restricted  to  be  source  free  so 
that  the  only  source  introduced  will  be  the  incident  plane  wave.  Therefore,  the  total  fields 
in  D  and  its  reflection  in  the  upper  half  plane  Di  satisfy  the  following  Maxwell’s  equations 
for  homogeneous,  linear,  isotropic,  source-free  media.  Using  the  time  convention,  these 
equations  are  given  as 


V  X  E  =  —jojfjiA 

(1) 

V  X  H  =  jweE 

(2) 

o 

II 

> 

(3) 

VH  =  0 

(4) 

where  (jj  is  the  radian  frequency  (radians/second),  fi  is  the  magnetic  permeability  of  the 
medium  (Henrys/meter),  e  is  the  electric  permittivity  of  the  medium  (Farads/meter),  and 
j  =  is  the  imaginary  unit  [6].  Taking  the  curl  of  equations  (1)  and  (2)  and  substituting 
in  appropriate  terms  reveals  that  both  E  and  H  satisfy  the  homogeneous  wave  equation 

V  X  V  X  E(f)  -  fc2E(r)  =  0  (5) 

V  X  V  X  H(f)  -  =  0  (6) 

where  k  =  is  the  propagation  constant  or  wave  number  of  the  medium  in  inverse 

meters. 
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In  addition  to  satisfying  Maxwell’s  equations  for  a  source-free  region,  the  total  fields 
must  also  satisfy  certain  boundary  conditions.  Since  the  ground  plane  and  the  smface  of 
the  cavity  S  are  perfect  electric  conductors,  the  tangential  component  of  the  total  electric 
field  will  be  “shorted  out”  there.  So  the  electric  field  must  satisfy  the  boimdary  conditions 


n  X  E(r)  =  0 

for  ?  G  5 

(7) 

z  X  E(f)  =  0 

for  r  G  cr*’ 

(8) 

where  is  the  complement  of  the  aperture  in  the  ground  plane  (i.e.  the  part  of  the  groimd 
plane  that  is  undisturbed  by  the  presence  of  the  cavity).  Furthermore,  the  tangential 
components  of  the  total  electric  field  and  the  total  magnetic  field  must  be  continuous 
across  the  aperture  a.  This  can  be  written  as 


lim 


z  X  E(r  -I-  Sz)  —  z  X  E(r  —  Sz) 


lim 


z  X  H(r -1-  (5z)  -  z  X  H(r  -  <5z) 


=  0 

for  r  G  <7 

(9) 

=  0 

for  r  G  <T 

(10) 

The  final  field  conditions  that  must  be  satisfied  pertain  to  the  scattered  fields.  Specifically, 
the  Sommerfeld  radiation  conditions  must  be  satisfied  at  infinity 


lim  r  •  VE®(r)  +  jA:|r|E®(r)l  =0  (11) 

|f|-foo  L  J 


lim  f  ■  VH®(?)  +  jfc|f  |H®(f)  1  =  0 

lr|-»-oo  L  J 


as  well  as  the  regularity  conditions  [26]: 


(12) 


lim  |r|E®(f)  =  0 

(13) 

|r|-»*oo 

lim  |r|H^(?)  =  0 

r  ->oo 

(14) 
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2-4  Green’s  Functions 


The  three-dimensional  scalar  Green’s  function  commonly  associated  with  electromag¬ 
netic  radiation  and  scattering  problems  is  defined  as 


G{k;T,r')  = 
R=\t- 


-jkR 


4irR 

f'l 


(15) 


where  r  =  xx  +  yy  +  zi  and  r'  =  x'x  +  y'y  -H  z'x  are  position  vectors  in  unprimed  and 
primed  coordinates,  respectively.  (This  will  be  converted  to  the  two-dimensional  case  and 
applied  to  the  cmrent  problem  later  in  the  thesis.)  In  [26],  the  Green’s  function  in  equation 
(15)  was  shown  to  satisfy  the  distributional  differential  equation 


V^G  +  k'^G  = -6{t  -  r') 


(16) 


where  5  is  the  Dirac  delta  distribution.  The  Green’s  function  also  satisfies  the  scalar 
analogues  of  the  radiation  and  regularity  conditions  given  in  the  previous  section. 

Define  the  auxiliary  scalar  functions 


Gp(fc;  f,  f)  =  G(A:;  F,  F')  -  G(A:;  F,  F^)  (17) 

GN{k-,  F,  F')  =  G(fe;  F,  F')  -h  G(fc;  F,  F])  (18) 

as  linear  combinations  of  the  scalar  Green’s  function  in  equation  (15)  with  different  argu¬ 
ments.  The  function  Gd  is  refered  to  as  the  Dirichlet  scalar  half-space  Green’s  function, 
and  Gn  as  the  Neumann  scalar  half-space  Green’s  function.  Next,  define  the  dyadic  func¬ 
tion 

F(fc;  F,  F')  =  -jkVG{k]  F,  F')  x  f  (19) 

where  I  =  xx  +  yy  -H  zz  is  called  the  idemfactor. 

Finally,  define  the  auxiliary  dyadic  functions 

Fi  =  -jk  (vGn  X  ft  -H  VGp  X  zz)  (20) 
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and 


r2  =  -jk  (vGd  X  It  +  VGtv  X  zzj  , 

(21) 

where  It  =  xx  +  yy  is  called  the  transverse  idemfactor.  It  can  be  shown  that  the  auxiliary 
dyadic  functions  in  equations  (20)  and  (21)  satisfy  the  homogeneous  boundary  conditions 

z  X  Fi  =  0  for  =  0 

(22) 

z  X  V  X  r2  =  0  for  z  =  0 

(23) 

as  well  as  the  equations 

z  •  [!(?)  X  V  X  Fi(fc;F,f')]  =  2z  •  [A(f)  x  V  x  ?(&;?,?')] 

for  z  =  0  (24) 

z  •  I^X(r)  X  r2(A:;r,f')]  =  2z  •  |^A(r)  x  r(A:;f,f')j 

for  z  =  0  (25) 

The  scalar  and  dyadic  functions  presented  in  this  section  were  first  introduced  in 
[4],  and  then  in  [26].  They  are  the  fundamental  building  blocks  in  the  development  of  the 
coupled  vector  integral  equations  for  the  scattering  problem  which  are  presented  in  the 
next  chapter. 
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IIL  Integral  Equation  Theory 

In  this  chapter,  the  fundamental  theorems  and  the  coupled  set  of  vector  integral  equations 
developed  by  Asvestas  and  Kleinman  in  [4]  and  extended  by  Wood  in  [26]  are  presented. 
The  set  of  vector  integral  equations  is  converted  to  a  set  of  scalar  integral  equations  for  the 
two-dimensional  scattering  problem  imder  TE  polarization  which  will  be  solved  approxi¬ 
mately  by  a  Method  of  Moments  algorithm. 


Fundamental  Theorems 

All  of  the  theorems  presented  in  this  section  first  appeared  in  [4]  and  were  later 
extended  in  [26].  The  proof  of  the  theorems  can  be  found  in  [26].  The  theorems  are  as 
follows. 

Theorem  1  Let  V  he  a  homogeneous  region  with  regular  boundary  dV,  n  be  the  outward 
unit  normal  vector  on  dV,  and  T  =  r(A;;f,  r^)  be  defined  as  in  equation  (19).  If  A  =  A(f) 
satisfies  V  x  V  x  A  —  k^A  =  0  V  r  eV,  then 


n  •  ^  X  (v  X  r) 


+  X  A  j  X  r  ds  =  < 


jkV  X  A(r') 
0 


t'eV 


r'  ^  F 


(26) 


where  V  is  the  closure  of  V;  that  is,  the  union  of  V  and  its  boundary. 

Theorem  2  Let  V,  V,  dV ,  n,  and  A  be  defined  as  in  Theorem  1.  Let  Fi  6e  defined  by 
equation  (20),  Vi  be  the  image  ofV  with  respect  to  the  xy-plane,  and  Vt  be  the  closure  of 
Vi.  Then 


J  n  •  A  X  X  Fi^  +  X  A^  X  Fi  ds  =  < 


;A;V'  X  A{r') 


r'€F 


jA:[v^x  A(r^)l.  f '  €  VJ 

J  I 


(27) 


t'  iVUVi 


11 


Theorem  3  Let  V,  V,  dV,  Vi,  Vi,  h,  ani  A  he  defined  as  in  Theorem  2.  Let  r2  be  defined 
by  equation  (21).  Then 


J  n-[Ax  (vxTa) +  (Vx  a) 


jkV  X  A(rO  r'eV 

-:/A:[v^x  A(r;)].  F' €  Fi  (28) 
0  r'^VUVi 


In  the  preceding  theorems,  the  del  operator  V  differentiates  with  respect  to  unprimed 
coordinates,  while  V'  differentiates  with  respect  to  primed  coordinates.  Also,  the  operator 
is  the  same  as  V'  except  that  ^  is  replaced  by  It  is  pointed  out  in  both  [4]  and 
[26]  that  the  theorems  do  not  hold  when  f '  lies  exactly  on  the  region’s  boundary  dV,  but 
that  they  are  valid  as  r '  approaches  dV  from  either  the  interior  or  the  exterior  of  V. 

The  theorems  relate  a  vector  function  evaluated  at  a  point  to  an  integral  over  a  closed 
surface.  Also,  notice  that  they  are  stated  in  very  general  terms.  The  only  requirements 
are  that  the  vector  function  A(r)  satisfy  the  homogeneous  wave  equation  and  that  the 
outward  pointing  rmit  normal  vector  n  be  defined  almost  everywhere  on  dV. 


3.2  Applications  of  the  Fundamental  Theorems 


The  theorems  presented  in  the  previous  section  can  now  be  applied  to  the  embedded 
cavity  scattering  problem.  This  is  accomplished  by  associating  V  with  either  the  cavity 
volmne  D  or  the  upper  half  space,  z  >  0,  and  by  letting  A(r)  take  on  the  various  field 
components  of  the  cavity  scattering  problem. 


Application  1  Let  F  =  D,  A  =  H,  and  r'  €  £>  in  Theorem  1.  By  applying 
Maxwell’s  equation  (2)  and  the  boundary  condition  in  equation  (7),  Theorem  1  can  be 
rewritten  as 

ifc?yiE(r')  =  - 

-jhYi  I  (zxE(r))-r(fci;?,rOdu. 


J  (n  X  H(r))  •  V  x  r(A:i;  r,  r')ds 
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In  equation  (29),  Yi  =  y/ei/m  is  the  intrinsic  admittance  of  the  cavity  material.  This 
application  shows  that  the  total  electric  field  at  any  point  in  the  cavity  volume  is  a  function 
of  the  tangential  electric  and  magnetic  fields  on  the  cavity  surface  and  aperture. 


Application  2  Let  F  =  D,  A  =  E,  and  r'  G  D  in  Theorem  1.  By  applying 

Maxwell’s  equation  (1)  and  the  bormdary  condition  in  equation  (7),  Theorem  1  can  be 

rewritten  as  _ 

k\Z{k{r')=  J  (zxE(f))-V  xr(fei; 


-jkiZi  ^  (n  X  H(r))  •  T{ki;?,f')ds 


(30) 


where  Zi  =  \fiiifei  is  the  intrinsic  impedance  of  the  cavity  material.  This  application 
demonstrates  that  the  total  magnetic  field  at  any  point  in  the  cavity  volume  is  a  function 
of  the  tangential  electric  and  magnetic  fields  on  the  cavity  surface  and  aperture. 

Application  3  Let  V  be  the  upper  half  space  («  >  0),  A  =  E®,  and  >  0  in 
Theorem  2.  By  applying  Maxwell’s  equation  (1)  and  the  conditions  in  equations  (8),  (22), 
and  (24),  Theorem  2  can  be  rewritten  as 


H®(? 


'>=/ 

•/O’ 


zxE(f)-Vxr(A;o;?,fOd<r. 


(31) 


According  to  the  third  application,  the  scattered  magnetic  field  at  any  point  in  the  upper 
half  space  is  a  function  of  the  tangential  electric  field  on  the  cavity  aperture. 

Application  4  Let  V  =  D,  A  =  E,  and  r'  G  D  in  Theorem  3.  Using  Maxwell’s 
equation  (1)  and  the  conditions  in  equations  (7),  (23),  and  (25),  Theorem  3  can  be  rewritten 
as 

yA:iH(r ')  =2  J  (z  x  H(r))  •  f  (fci;  r,  r')da 
+  J  ^n  X  H(f)^  •r2(fci;r,f')ds 

This  application  shows  that  the  total  magnetic  field  at  any  point  in  the  cavity  voliune  is  a 
function  of  the  tangential  magnetic  field  on  the  cavity  surface  and  aperture. 

Application  5  Let  V  be  the  upper  half  space  (z  >  0),  A  =  H®,  and  z'  >  0  in 
Theorem  3.  By  applying  Maxwell’s  equation  (2)  and  the  conditions  in  equations  (8),  (23), 


(32) 
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and  (25),  Theorem  3  can  be  rewritten  as 

— ^E®(r ')  =  f  zx  E(f)  •  r(fco;  r,  r')da.  (33) 

2  Jcr 

According  to  this  application,  the  scattered  electric  field  at  any  point  in  the  upper  half 
space  is  a  function  of  the  tangential  electric  field  on  the  cavity  aperture.  Equation  (33) 
will  be  used  later  to  determine  the  monostatic  RCS  values  once  the  imknown  tangential 
field  components  have  been  found. 

The  applications  of  the  fundamental  theorems  define  the  total  fields  in  the  cavity 
volume  and  the  scattered  fields  in  the  upper  half  space  in  terms  of  the  tangential  field 
components  on  the  cavity  surface  S  and  the  cavity  aperture  a.  So  it  is  sufficient  to 
find  these  tangential  field  components  to  solve  the  cavity  scattering  problem.  Examining 
equations  (29)  through  (33),  it  is  apparent  that  there  are  two  different  types  of  field 
quantities  involved.  The  left  hand  sides  of  equations  (29)-(33)  involve  field  quantities 
evaluated  at  points  away  from  the  boundary  of  the  cavity  dD^  whereas  the  right  hand 
sides  of  the  equations  involve  the  tangential  components  of  the  total  fields  evaluated  on 
dD. 

In  order  to  make  equations  (29)“(32)  more  useful,  consider  letting  r'  approach  dD 
and  then  sifting  out  the  tangential  component  of  the  result  as  suggested  by  Wood  in  [26]. 
Begin  by  fixing  a  point  f  on  dD  so  that  n  is  the  outward  unit  normal  vector  at  r,  and 
then  cross  multiply  the  equations  by  n.  Finally,  let  r '  —  r  +  5n  and  evaluate  the  limit  as 
(5-^0.  Wood  [26]  also  presents  a  theorem  to  help  facilitate  the  limiting  process. 

Theorem  4  Let  V  be  a  volume  with  regular  boundary  dV ,  f  be  a  fixed  point  on  dV ,  n  be 
the  outward  unit  normal  vector  atv,  and  r(A;;r,  r')  be  defined  as  in  equation  (19).  If  A{v) 
is  continuous  on  dV,  then 

limnx  [  l(f")-r(A:;r",f')rfs"  =  T^-i(?)  +  Ax  [  A{f")  -  (34) 

r'-4r  Jqv  JaV 

where  the  upper  sign  is  taken  ifr'—^r  from  the  exterior  of  V,  and  the  lower  sign  is  taken 
iff'  f  from  the  interior  ofV. 
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Application  la  Fix  r  G  (t  so  that  n  =  z,  and  set  r'  =  r  —  (5z,  where  (5  >  0.  Cross 
multiplying  equation  (29)  by  z  and  taking  the  limit  as  5  0  via  Theorem  4  simplifies  the 

equation  to 

^^^^z  X  E(r')  =zx  f  fn  X  H{r))  •  V  x  v')ds  for  f '  G  cr.  (35) 

2  JdD  '  ^ 


Notice  that  the  second  integral  in  equation  (29)  vanished  in  the  limiting  process.  This  is 


because  z  x  ^z  x  E(r)^  r,  ?') 


=  0  when  f  G  a  and  f '  G  <t.  This  new  application 
shows  that  the  tangential  electric  field  on  the  aperture  a  can  be  expressed  as  a  function  of 
the  tangential  magnetic  field  on  dD. 


Application  2a  Fix  r  G  <r  so  that  n  =  z,  and  set  r'  —  r-  Sz,  where  <5  >  0.  Again, 
cross  multiplying  equation  (30)  by  z  and  taking  the  limit  as  (5  — 0  using  Theorem  4,  the 
equation  can  be  written  as 


^l^z  X  H(r')  =  z  X  J  •  V  X  T{ki;r,r')da 

—jkiZiz  X  J  ^nx  H(r)^  •  r(A:i;r,f')ds  for  r'  G  <t. 


(36) 


This  application  now  says  that  the  tangential  magnetic  field  on  a  can  be  expressed  as  a 
function  of  the  tangential  magnetic  field  on  S  and  the  tangential  electric  field  on  a . 

Application  3a  Fix  r  G  cr  so  that  n  =  z,  and  set  r'  =  r  +  6z,  where  (J  >  0.  Cross 
multiply  equation  (31)  by  z  and  take  the  limit  as  <5  — >•  0  using  Theorem  4  to  obtain 


— ^^^z  X  H®(r')  =  z  X  fix  E(f)  •  V  x  r(feo;  f,  r')da  for  r '  G  cr.  (37) 
2  Jo’ 

Recall  that  H®  =  H— Also  note  that  along  the  aperture  since  z  xH*“®  =  z  x 
equation  (37)  can  be  rewritten  as 


z  X 


[H(f ')  -  2H*"®(f ')]  = 

X  r(fco;r,f')cicr  for  r' G  cr. 


(38) 


15 


Application  4a  Fix  r  e  5  and  set  r'  =  r-5z,  where  (5  >  0.  Cross  multiply  equation 
(32)  by  n'  and  take  the  limit  as  5  — 0  using  Theorem  4  to  obtain 


^n'  X  H(rO  =  2n'  X  j  (zx  H(f))  •  r(fci;r,r')dcr 
+n'  X  J  (nx  H(r)j  •  r2(fci;r,r')ds  for  r'  E  S. 


(39) 


The  equations  in  Applications  2a  and  3a  can  be  combined  and  simplified.  Doing  so 
results  in  a  system  of  three  coupled  vector  integral  equations  that  govern  the  electromag¬ 
netic  scattering  from  an  embedded  cavity  in  a  perfectly  conducting  ground  plane.  Formally, 
the  system  of  vector  integral  equations  can  be  written  as 


X 


y  X  H(r) j  •  V  X  r(fci; f, r')d<T  -h  z  x  J  x  H(f) j  •  V  x  T{ki;r, ?')ds 


=  X  E(rO  for  F'  €  a, 


(40) 


d<T  X  H(r  0 


-Ziz  X  [  hx  H(r)  •  ?,  ?')ds  =  jIzqZqz  x  H'“"(r ')  for  f '  G  a, 

Js 


(41) 


and 


2n'  X  J  (zx  H(r))  •  r(A:i; r, r')d<T  -I-  n'  x  J  (n  x  H(r))  •  r2(A:i; r, r')ds 


=  ^n'  X  H(r')  for  r'  G  S. 


(42) 


3.3  Derivation  of  Scalar  Integral  Equations 

In  this  section,  the  vector  integral  equations  (40)- (42)  are  converted  to  a  system  of 
scalar  integral  equations  for  the  scattering  problem,  where  the  incident  plane  wave  is  under 
transverse  electric  (TE)  polarization.  The  derivation  is  similar  to  the  method  presented 
in  [27]  where  scalar  equations  were  developed  for  transverse  magnetic  (TM)  polarization. 
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The  necessary  field  quantities  and  unknowns  are  defined  and  then  substituted  into  the 
vector  equations.  These  are  then  simplified  and  converted  to  scalar  equations  which  can 
be  solved  approximately  by  a  Method  of  Moments  (MoM)  algorithm. 

Under  TE  polarization,  all  of  the  components  of  the  magnetic  field  are  parallel  to 
the  longitudinal  x-axis.  For  a  plane  wave,  define  the  magnetic  portion  of  the  incident  field 
to  be 

z)  =  e^ko(!/cose+zsinO)^  (43) 

and  the  unknown  total  magnetic  field  as 

H(y,z)  ^u(y,z)x  (44) 

where  u(y,z)  is  some  unknown  continuous  function.  Then  applying  Maxwell’s  equations, 

the  corresponding  incident  and  total  electric  fields  can  be  shown  to  be 

z)  =  Zo(sin0y  -  (45) 

and 

E(?/,  z)  = - Vu  X  X.  (46) 

jcoei 

It  should  be  pointed  out  here  that  the  vector  integral  equations  in  the  previous  section 
are  valid  in  three  dimensions  as  well  as  in  the  two-dimensional  problem  now  considered. 
The  only  diflFerences  are  that  the  field  quantities  are  no  longer  functions  of  the  longitudinal 
variable  x  and  that  the  appropriate  two-dimensional  Green’s  function  must  be  used  instead 
of  the  three-dimensional  Green’s  function.  This  will  be  defined  when  the  scalar  equations 
are  solved  in  a  later  section.  For  notational  purposes,  the  two-dimensional  position  vector 
p  =  yy  +  zi  will  now  be  used  in  regard  to  the  field  quantities.  Doing  so,  the  tangential 
field  components  appearing  in  equations  (40)-(42)  can  be  defined  as  follows. 

z  X  H(p)  =  u{^y  for  p  €  cr  (47) 

n  X  H(^  =  u(p)n  x  x  for  p  G  S'  (48) 
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z  X  E(/o)  = - —  for  ^  G  <T  (49) 

juei  az 

=  forpGtr  (50) 

Now  these  quantities  can  be  substituted  into  the  vector  integral  equations,  whidi  can  then 
be  simplified  to  produce  scalar  integral  equations. 

3.3.1  Scalar  Equation  1.  The  only  other  term  appearing  in  equation  (40)  that 
still  needs  to  be  defined  is  the  ciu’l  of  the  dyadic  function  F.  Recall  that  it  was  of  the  form 

F  =  -jkVG  X  f  (51) 


For  readability,  the  arguments  of  the  Green’s  functions  and  of  the  dyadic  functions  will 
be  suppressed  for  the  remainder  of  this  chapter.  Also,  in  order  to  indicate  which  wave 
nmnber  appears  in  the  argument  of  the  Green’s  functions,  it  will  be  understood  that 
Gi  =  G{ki]  p,  p')  and  Go  =  G(fco;  p,  p').  Given  this  notation,  it  can  be  shown  that  the  cm:l 
of  the  dyadic  function  in  equation  (40)  can  be  written  as 

V  X  F  =  -]ki  (kfGii  +  VVGi)  (52) 


where 


a^Gi^  d^Gi^  52Gi_  a^Gi. 

vvGi  = 


(53) 


dy2  '  9^2  dy9z 

is  the  Hessian  of  the  two-dimensional  Green’s  function.  Substituting  each  of  the  newly 
defined  components  into  the  first  vector  integral  equation  yields 


z  X  /  u{p)y  ■  {-jkx)  ffciGiI  +  VVGi j  da 

J  (T 

+z  X  ^  («(/^n  X  x)  •  i-jki)  (kfGj.  +  VVGi)  dl  = 


(54) 


Denote  the  components  of  the  newly  formed  integral  equation  as  I  +  II  =  III.  To  begin 
simplifying  the  first  part,  take  the  dot  product  with  y  as  follows. 
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Now  taking  the  cross  product  of  z  with  the  integrand  reduces  I  to 


/  =  X 


L 


3klu{^Gi+  jkiu{^ 


dy^  J 


da. 


(56) 


The  second  part  of  the  integral  equation  can  be  written  as 


77  =  z  X  J  —jkiu{^h  •  X  X  +  VVGi^j  dl 


(57) 


where  the  identity  (a  x  b)  •  c  =  a  •  (b  x  c)  has  been  used  [9].  Next,  performing  the  cross 
product  within  the  integrand  yields 


77  =  z  X  /  — ^fci«(p)n  • 

Js 


a"G,..  a=Gi..  ,  ,  a^'Ci.: 

*,G,(^y  -  y.)  -  -^yz  -  ^yy  +  -^zy  +  ^zz 


dl. 

■(58) 


Dotting  n  into  the  integrand  and  then  taking  the  cross  product  of  z  with  the  integrand 
reduces  it  to 


II  =  x  J^jkiui^  nz(klGi  +  ^^^ 


dl 


(59) 


Qy2  J  '"^dydz 

where  Uy  and  %  are  the  appropriate  components  of  the  outward  pointing  imit  normal 
vector  along  S.  Finally,  the  right  hand  side  of  the  first  vector  integral  equation  can  be 
written  as 

(60) 


‘2j  dz 

Dotting  X  into  both  sides  of  the  integral  equation  represented  by  the  newly  formed  I + II  = 
III,  dividing  both  sides  by  jki,  and  simplifying  reduces  equation  (40)  to 


j  u{^  (k^Gi  +  is 


2^  . 


Gi  + 


%2 


dydz\ 


dl 


(61) 


Equation  (61)  is  now  a  scalar  integral  equation  which  is  in  a  form  that  can  be  solved 
approximately  by  the  numerical  MoM  algorithm. 
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3.3.2  Scalar  Equation  2.  Substituting  the  field  components  into  the  second  vector 
integral  equation  resiilts  in 

-z  X  /  [(fc?Gif+ VVGi)  -  (A:gGoi  + VVGo)]  da 

-ZiixJ^  {u{^h  X  x)  •  (-jfeiVGi  X  i)  dl  (62) 
=  for  p'  G  <t. 

As  before,  break  equation  (62)  into  fomr  parts  so  that  it  is  of  the  form  I +  II  +  III  =  IV. 
Then  the  first  part  can  be  written  as 

J  =  z  X  f  •  [f  [klGi  -  klGo)  +  VV  (Gi  -  Go)]  da.  (63) 

dz  i  ^  J 

Taking  the  dot  product  within  the  integrand  reduces  the  equation  to 

/  =  zx  f  ^  (fcf  Gi  —  fcpGo)  X  da.  (64) 

jujey  dz  ^ 

Notice  that  x  •  [VV  (Gi  —  Go)]  vanished  in  the  previous  step.  This  is  due  to  the  fact  that 
the  Green’s  functions  Gi  and  Go  have  no  x  components,  so  consequently,  the  Hessian  of 
Gi  —  Go  has  no  x  component.  Finally,  taking  the  cross  product  of  z  with  x  in  7  reduces 
it  to 


Part  7/  of  equation  (62)  is  already  in  simplified  form.  The  third  part  can  be  reduced  by 
first  performing  the  cross  multiplication  of  the  gradient  of  Gi  with  the  idemfactor  to  obtain 

777  =  jkiZiz  X  J  («(p)n  x  x)  •  ~ 

Next,  taking  the  dot  product  of  x  with  the  right  side  of  the  integrand,  followed  by  cross 
multiplying  the  integrand  with  n  simplifies  part  777  to 

III  =  jkiZiix 
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The  kernel  of  equation  (67)  is  easily  recognized  as  the  normal  derivative  of  Gi-  So,  cross 
multiplying  z  with  the  integral  simplifies  it  to 


III  =  yjkiZi 


(68) 


The  right  hand  side  of  equation  (62)  is  already  in  simplified  form,  so  the  final  step  is  to 
take  the  dot  product  of  y  with  both  sides  of  the  equation  represented  by  /+// +III  =  IV. 
Doing  so  gives  the  second  simplified  scalar  integral  equation  as 


J_  r  {klGi-klGo)da  +  jhZi  f 
juei  Ja  dz  ^  Js  on 

+fw(pO  (feiZi  +  koZo)  =  p'  6  tr. 


3.3.3  Scalar  Equation  3.  Substituting  the  necessary  components  into  equation 
(42)  simplifies  it  to 


2n'  X  J  «(p)y  •  jfeiVGi  x  da+ 


h'x  f  (u(^n  X  x)  •  [{-jki)  (vGd  X  It  +  VGn  x  zz)]  dl  =  ^u(p')n'  x  x,  p '  G  5 
^  (70) 

where  It  =  jcx+yy,  and  Go  and  Gn  are  defined  as  in  equations  (17)  and  (18),  respectively. 
In  a  similar  fashion  as  before,  the  first  part  of  the  equation,  denoted  as  7,  can  be  rewritten 


as 


7  =  2h' 


[  ui^y- 

J  <X 


i-jki) 


xz 


dGi 

dz 


xy  + 


- — zx  I  da. 

dy  )\ 


(71) 


Performing  the  dot  multiplication  then  gives 


7  =  -2jfci(n'  X  x) y  u{^^^da 


(72) 


as  the  simphfied  version  of  7.  The  second  part  of  the  integral  equation  can  be  written  as 


77  =  —jkixJ  X  [  {u{^n  x  x) 

Js 


-dGo..  ,  dGn..  BGb..  , 

— xy  +  — yx  -  — zx  +  — xz 


dl  (73) 
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where  the  required  cross  products  have  been  performed.  Again,  applying  the  identity 
(a  X  b)  •  c  =  a  •  (b  X  c),  the  second  part  of  the  integral  equation  can  be  reduced  to 

II  =  -jkivlx  J  u(p)n-  ^^yx+^^zx  dl.  (74) 

Recognizing  the  kernel  of  the  integrand  in  equation  (74)  as  the  gradient  of  Gp  being  dotted 
with  n  allows  the  second  part  to  be  simplified  as 

II  =  -jfci(n'  X  x)  y*  (75) 

Since  the  third  part  of  the  third  integral  equation  is  aheady  simplified,  the  only  thing 
left  to  do  is  to  dot  the  vector  (n'  x  x)  into  both  sides  of  the  equation  represented  by 
/  +  //  =  III  and  divide  both  sides  by  —jki.  This  gives  the  final  simplified  form  of  the 
third  scalar  integral  equation  as 

2  ju{p)^dl  =  ^u{p\  p'es.  (76) 

Equations  (61),  (69),  and  (76)  represent  the  scalar  form  of  the  three  coupled  vector 
integral  equations  that  govern  the  electromagnetic  scattering  from  an  embedded  cavity  in 
a  PEC  ground  plane  imder  TE  polarization.  The  integral  equations  are  now  in  a  form 
which  can  be  approximately  solved  for  the  unknown  fields  using  a  numerical  technique. 
The  numerical  MoM  algorithm  is  developed  in  the  following  chapter. 
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IV.  Numerical  Solution 


In  this  chapter,  the  general  procedure  of  solving  a  linear  operator  equation  via  the  Method 
of  Moments  (MoM)  is  discussed.  This  is  followed  by  the  development  of  a  MoM  algorithm 
designed  to  find  an  approximate  solution  to  the  three  scalar  integral  equations  (61),  (69), 
and  (76).  In  applying  the  MoM  algorithm,  certain  singular  integrals  arise  and  must  be 
handled  properly.  This  chapter  demonstrates  how  such  integrals  can  be  treated  so  that 
they  can  be  evaluated  numerically.  Finally,  the  error  of  the  MoM  algorithm  is  considered 
and  an  estimate  of  the  order  of  convergence  of  the  method  is  given. 

4-1  Method  of  Moments 

The  basic  premise  behind  the  Method  of  Moments  is  to  transform  a  linear  operator 
equation  into  a  system  of  simultaneous  linear  equations.  The  method  was  first  conceived 
by  a  Russian  scientist  named  Galerkin  in  1915  to  solve  elastodynamic  problems.  Since 
then  many  scientific  disciplines  have  utilized  the  idea  to  find  an  approximate  solution 
to  problems  involving  linear  integral  and  differential  operators.  While  many  scientists 
contributed  to  the  development  of  the  theory,  in  1967  Roger  Harrington  of  Syracuse  Uni¬ 
versity  presented  a  unified  approach  to  solving  electrod3mamic  problems,  including  the 
electromagnetic  scattering  problem,  in  a  generalized  moment  method  formulation  [8]. 

The  MoM  is  applied  to  problems  in  which  the  unknown  quantity  is  embedded  in  an 
integrand  or  is  the  argument  of  a  differential  operator  and  cannot  be  solved  for  analytically. 
The  process  of  a  linear  operator  acting  upon  the  unknown  to  produce  a  functional  effect 
can  be  written  mathematically  as 

L[u]  =  f,  (77) 

where  L  is  the  linear  operator,  u  is  the  unknown  function,  and  /  is  the  effect  or  the  forcing 
function.  The  unknown  quantity  u  is  first  expanded  in  terms  of  a  chosen  basis  function 
set,  that  is 

00 

Y.  (78) 

n,=— oo 
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where  the  In  are  unknown  expansion  coefficients  and  the  axe  functions  from  the  chosen 
basis  set.  Obviously,  an  infinite  number  of  terms  cannot  be  used  in  the  computation  of 
an  approximate  solution  to  equation  (77),  so  the  series  is  trimcated  to  a  finite  number  of 
terms  plus  a  final  term  which  represents  all  of  the  truncated  terms,  i.e.  the  error  e.  This 
is  shown  by  the  equation 

N 

u  =  ^  Inipn  +  e.  (79) 

Substituting  the  above  expression  for  u  into  equation  (77)  gives 


L 


-  N 


Ln=l 


=  /, 


which  due  to  the  linearity  of  L  can  be  written  as 


(80) 


N 

In  {L  [ipn])  +re  =  f- 

n=l 


(81) 


In  equation  (81),  rg  is  called  the  residual.  For  an  accurate  approximation  it  is  desirable  to 
minimize  the  residual.  This  is  accomplished  by  taking  an  inner  product  with  a  family  of 
testing  functions,  0^)  which  are  forced  to  be  orthogonal  to  the  residual.  That  is. 


(®m)  —  0, 


(82) 


where  the  inner  product  is  defined  to  be 


{em,A.)  =  Jq 


0m  •  A  dl. 


(83) 


To  perform  the  testing  procedure,  the  inner  products  of  both  sides  of  equation  (81)  are 
taken  with  the  testing  functions  which  gives 


N 

^  ]  In  [^n] )  ®m)  d"  —  (®t7u  /)  • 

n=l 


(84) 
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Applying  the  orthogonality  condition  of  equation  (82)  and  rearranging  terms  results  in 

N 

{em,f)  =  J2ln{L[^n],em),  for  m  =  1, 2, . . .  , Af.  (85) 

n=l 

For  the  choice  of  N  basis  functions  and  N  testing  functions,  the  equations  in  (85)  represent 
a  system  of  N  linear  equations  with  N  unknowns,  the  expansion  coeflBcients  In-  This  can 
be  written  in  matrix  notation  as 


[V]n>.i 

—  [^]nxN  [^liVxl  ’ 

(86) 

where 

=  (0m,/) 

[^]mn 

=  (i[V’n],0m) 

(87) 

mn 

=  In 

for 

m 

n  ■ 

II  II 

(88) 

In  equation  (87),  [V]^  is  often  called  the  voltage  vector  or  the  excitation  vector,  [Z]^^  is 
referred  to  as  the  impedance  matrix  or  the  coupling  matrix,  and  [7]^  is  called  the  cmrrent 
vector  or  the  solution  vector  [10]. 

Oftentimes  the  integrations  involved  in  calculating  the  elements  of  the  coupling  ma¬ 
trix  are  difficult  to  perform,  so  it  is  desirable  to  choose  a  set  of  testing  functions  that  will 
simplify  the  process.  One  of  the  most  commonly  used  techniques  to  obtain  approximate 
solutions  to  these  integrations  is  called  point  matching.  The  point  matching  procedure 
simply  requires  that  equation  (81)  be  satisfied  exactly  at  discrete  points,  called  match- 
points,  in  the  region  of  interest.  This  is  accomplished  by  using  Dirac  delta  functions  as 
testing  functions.  In  doing  so,  the  inner  product  in  equation  (83)  boils  down  to  a  function 
evaluation  which  can  be  easily  performed  [12]. 
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4-2  Cavity  Discretization 


In  order  to  find  approximate  solutions  to  the  scalar  integral  equations  (61),  (69),  and 
(76),  the  cavity  perimeter  dD  —  a  U  S  is  discretized  into  non-overlapping  segments.  Let 
there  be  a  total  of  Ni  segments  along  the  cavity  aperture  a  and  a  total  of  N2  segments 
along  the  cavity  surface  S.  The  segments  are  defined  by  the  node  points  {yk^iUkih)  as  ™ 
Figure  3.  For  each  node,  {yki^k)  represent  the  (y,  z)  coordinates  of  the  point,  and  Ik  is  the 
arc  length  measured  clockwise  from  (yi,  zi)  to  (y*,  Zk)  along  dD.  In  the  example  shown  in 
Figure  3,  IVi  =  2  and  N2  =  5. 


4.3  Method  of  Moments  Algorithm 

In  this  section,  the  MoM  is  applied  to  the  three  coupled  scalar  integral  equations 
(61),  (69),  and  (76)  to  find  approximate  solutions.  The  procedure  is  very  similetr  for  all 
three  equations,  so  some  of  the  details  have  been  omitted  from  the  development  of  the 
MoM  algorithm  for  the  second  and  third  equations  for  brevity.  Later  in  the  chapter,  an 
explanation  of  how  to  program  the  MoM  algorithm  on  a  computer  as  well  as  how  to  deal 
with  certeiin  singular  integrals  that  arise  is  given. 
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4-3.1  Solution  Scheme  for  Equation  1.  In  this  section  the  MoM  algorithm  is 
developed  for  the  first  scalar  integral  equation.  Recall  that  it  is  of  the  form 


J^u{^  (kfGi  +  da  +  - 


cf^Gi 


dydzj 


dl 


-1  du{p')  , 

p  €  a 


(89) 


dz 


where  Gi  is  the  appropriate  two-dimensional  Green’s  function.  The  unknown  quantities 
in  equation  (89)  are 

ui^  pea  (90) 


u{p') 

du{^ 

dz 


^eS 

pe  a 


(91) 

(92) 


which  represent  the  total  magnetic  field  along  the  cavity  aperture,  the  total  magnetic  field 
along  the  cavity  surface,  eind  the  normal  derivative  of  the  magnetic  field  along  the  cavity 
aperture,  respectively. 

The  first  step  in  the  MoM  procedure  is  to  expand  the  unknowns  in  terms  of  the 
chosen  basis  function  set.  Using  pulse  basis  functions,  the  unknowns  in  equations  (90), 
(91),  and  (92)  can  be  approximated  as 


N1+N2 

u(0  ~  ^  ]  0,nPn{^)i  p{0  ^ 


n=l 


and 


du{l) 

dz 


Ni 

n—l 


(93) 


(94) 


where  a„  and  are  the  unknown  expansion  coefiicients  and  p„(/)  is  the  pulse  basis  ftmc- 
tion.  Note  that  the  vector  position  functions  are  now  defined  in  terms  of  the  arc  length  I 
along  the  perimeter  of  the  cavity.  The  pulse  basis  function  is  then  defined  to  be 


Pn{l) 


1 

0 


otherwise 


(95) 
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As  mentioned  earlier,  Dirac  delta  functions  will  be  used  as  the  testing  functions  to  employ 
the  point  matching  technique  and  simplify  calculations.  This  is  defined  to  be 


(96) 


where 


“1“  ^m+1 


(97) 


is  the  midpoint  or  matchpoint  of  each  segment. 

Now  the  approximations  in  equations  (93)  and  (94)  can  be  substituted  for  the  un¬ 
knowns  in  Equation  (89)  to  produce 


+ 


n=Ni+l 


Ni 


Tl=l  ^ 

Next  the  inner  product  used  will  be  defined  as 

rr 


(e„,A)  =  ^ 


(99) 


where  T  is  the  total  arc  length  around  the  perimeter  of  the  cavity.  Taking  the  inner  product 
with  &fn  of  both  sides  of  equation  (98)  gives 

E  »»  £  ['*'  pM  !')  +  ■«  (*'  -  *m+i)  «'+ 


Pn(0  1^"^  (/,/)+  Qy2  j  %  QyQ^  L+^)dldl 


Using  the  sifting  property  of  the  Dirac  delta  function,  equation  (100)  simplifies  to 


This  can  be  written  more  concisely  as 

Ni  iVi+N2 

I 

n=l  n=iVi+l 

where 

rln+l 


Ni  N1+N2  ^  Ni 

^  ^  H"  ^  ^  ^nPmn  d"  «  ^  ^  “  0 


n=l 


9S/2 


and 


Pmn  —  I  fe^Gi  +  ^^2 


- 


dydz 


for  m  =  1, 2, 3, ... ,  JVi  and  n  =  1, 2, 3, . . . ,  JVi  +  JV2.  Equations  (103)  and  (104) 
be  used  to  fill  the  entries  of  a  matrix  of  the  form 


^1  = 

I 

OiNi,l 

The  linear  system 


ai,i  •••  ai,JVi  ^i,Ni+N2  k  ® 

OiNuNi  I3ni,Ni+1  •••  Pni,Ni+N2  0  •••  ^ 


r 

1 

0 

_ 1 

■ 

1 

- 1 

•  0 

_ 1 

(101) 

(102) 

(103) 

I  (104) 

can  now 

(105) 

(106) 
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contains  2Ni  +  unknown  expansion  coefficients  but  only  Ni  equations.  The  remaining 
equations  necessary  to  obtain  a  unique  solution  will  come  from  the  other  two  scalar  integral 
equations. 


4.3.2  Solution  Scheme  for  Equation  2.  In  this  section  the  MoM  algorithm  for 
the  second  scalar  integral  equation  is  developed.  The  equation  is 


(fcfCi  -  klGo)  da  +  jkiZi  f 
juiei  Ja  oz  ^  Js  on 

+  ^~u{p')  (fciZi  +  koZo)  =  p'  6  (t. 


(107) 


Proceeding  in  a  similar  fashion  as  the  previous  section,  pulse  basis  functions  are  used 
to  expand  the  unknowns  in  terms  of  the  expansion  coefficients  an  and  hn-  Discretizing 
equation  (107),  taking  the  inner  product  with  the  delta  testing  functions,  and  rearranging 
terms  gives 


(i,UL)  -  *«G,  di 

^4+^2  A+i 

(fciZi  +  fco.^0)  +  ^  ~ 


dl 


Equation  (108)  can  be  written  as 

Ni  N1+N2  Ni 

H  On  +  0,nVmn  +  hn'tmn  =  /  ’ 

n=l  ri=iVi+l  n=l 


(108) 


(109) 


^  =  |(fci  ^1  +  4^0), 


Vmn 


=  3kiZi  [ 

Jin 


ln+1  dGl 

dn 


dl, 


7mn  —  ^^^1  ^  ^0^0  (^’ ^m+^)] 


(110) 

(111) 

(112) 
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where 


and 


/  =jhZoe 


(113) 


jkoiy, 


'”»+5 


COS  9-\~z 


^^,sin0) 


form  =  l,2,3,...,iVi  andn  =  1,2,3, ...  ,iVi +iV2.  Now  equations  (110),  (111),  and  (112) 
are  used  to  compute  the  elements  of  a  matrix  of  the  form 


A2  = 


H  •  •  0 


0  •••  H  77iVi,iVi41 


'n\,Ni+Ni  71,1  *  •  •  ll,Ni 


IJiVulVi+ATa  7lVi,l  •••  iNuNi 


(114) 


The  system  represented  by  A2  provides  an  additional  Ni  equations  with  the  same  unknowns 
as  the  system  in  (106).  The  third  integral  equation  will  provide  the  other  necessary  N2 
equations. 


4.3.3  Solution  Scheme  for  Equation  3.  In  this  section  the  MoM  algorithm  is 
developed  for  the  final  scalar  integral  equation  which  is  given  by  the  formula 

+  p'€S.  (116) 


Following  the  same  procedure  as  before,  equation  (115)  is  discretized  using  pulse  basis 
functions,  and  the  inner  product  is  taken  with  the  delta  testing  functions.  This  results  in 
the  equation 


A  /•*" 

n=l 


8G.(i,Uj)  "■-I"’ 

—dl  + 


dz 


4-1  •'‘n 


ln+1  dGoihljn+h) 


-1 


n=Ni+l 
ATi+JVa 

n=iVi+l 


dn 


dl 


(116) 


This  can  be  more  simply  written  as 

Ni  iVi+iVs  N1+N2 

^  ^  ^Ti4^mn  ~l~  ^  ^  "b  2  ^  ^  On  ”  O9  (117) 

Ti=l  n=Ni-\-l  n=Ni-\-l 
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where 


(118) 


and 


Jin 


Jn+1  dGl{l,lj^,l) 


- 

Jin 


dz 

In+i  dGoilJjn^i) 


dn 


-dl 


dl 


(119) 


for  m  =  iVi  +  1,  iVi  +  2, . . . ,  iVi  +  IV2  and  n  =  1, 2, 3, . . . ,  iVi  +  N2.  Using  these  equations, 
the  following  matrix  can  be  computed 


-As  = 


V’JVi-fl,!  •••  'lpNi+l,Ni  ^Ni+l,Ni+l  +  ^ 

'll’Ni+N2,l  •••  ‘fpNi+N2,Ni  ^Ni+N2,Ni+1 


^Ni+1,Ni+N2  0 


0 


•  •  ^Ni+N2,Ni+N2  +5  0  •  •  •  0  J 

(120) 

Notice  that  there  is  an  JV2  by  Ni  submatrix  of  zeros  appearing  on  the  right  side  of  As. 
This  is  because  the  unknown  in  (92)  does  not  appear  in  equation  (115).  Using  this  MoM 
algorithm,  the  coupled  set  of  integral  equations  has  been  converted  to  the  linear  system 
Au  =  /,  where 

Ai  A2  A^  j  5  (121) 

=  [  ai  ...  aiVi+iV2  h  ...  6isri  ]  >  (122) 


and 


/=[o  ...  0  /(;„p  ...  /(i„.+i)  0  ...  o]  . 


(123) 


4-4  Computing  the  Impedance  Matrix 

Now  that  the  MoM  algorithm  has  been  developed,  it  is  necessary  to  simplify  the 
equations  used  to  compute  the  elements  of  the  impedance  matrix  A.  For  example,  equa¬ 
tions  (103)  and  (104)  are  not  currently  in  a  form  that  can  be  easily  implemented  in  a 
computer  code  to  produce  the  elements  of  Ai.  By  applying  some  common  identities  and 
defining  some  of  the  expressions  in  the  MoM  algorithm,  the  elements  of  A  can  be  easily 
calculated. 
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Recall  the  three-dimensional  scalar  Green’s  function  was  defined  earlier  as 


g-jkR 

47ri2 


R=|r-f'| 


(124) 


where  r  =  x^  +  yy  +  zz  and  r'  =  x'x  +  y'y  +  2r'z  are  position  vectors  in  unprimed  and 
primed  coordinates,  respectively.  The  two-dimensional  Green’s  function  can  be  found  by 
making  use  of  the  identity 

/oo  ^—jkR  1  ,  . 

^dx  =  =  G‘^^{k;p,p'),  (125) 

where  is  the  Hankel  function  of  the  second  kind,  and  again,  p  =  r  •  (yy  -t-  zz)  is  the 
projection  of  the  position  vector  onto  the  j/2r-plane. 

For  each  segment  on  dD  defined  by  [/„,  /n+i]?  {Vni  ^n)  represents  the  first  node  point 
and  (yn+i^Zn+i)  represents  the  second  node  point.  Using  these  node  points,  the  compo¬ 
nents  of  the  position  vector  along  dD  can  be  defined  in  terms  of  the  arc  length  as 

y(0  =  yn  + ,  ^  S  (yn-t-1  -  yn)  (126) 


and 


41)  =Zn+  /  S  (Zn+l  -  Zn)  -  (127) 

*n+l  “ 

The  components  of  the  discrete  set  of  matchpoints,  which  are  the  midpoints  of  the  seg¬ 
ments,  can  be  expressed  as 


Vm  +  ym-f  1 
^"*+5  “  2 

^m+l  -  2  ■ 

Using  these  definitions,  the  difference  \p  —  p'\  in  equation  (125)  can  be  expressed  as 


(128) 
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Then  the  two-dimensional  Green’s  function  for  the  MoM  algorithm  can  be  expressed  as 


G  (fc;  I,  1 )  =  ^  [Jo  (k  •  dm{l))  -  jYo  {k  ■  dmil))] . 


(130) 


4.4.1  Computing  Ai.  Recall  the  amn  elements  of  Ai  are  given  by  the  equation 


CXr] 


t'  ‘icd/.Ui) 


+ 


dy^ 


dl. 


(131) 


This  can  be  rewritten  as 


Oimn 


)  (kidmil))  + 


1  {kidmjl)) 

Aj  dy^ 


dl. 


(132) 


Using  the  properties  of  Bessel  functions,  the  second  derivative  term  in  equation  (132)  can 
be  simplified  to 


(kidrajl))  _ 

dj/2 


kf 


-iff  (kidmil))  + 


1 


kidm{l)^^ 


(2) 


(kidmil)) 


(133) 


Plugging  this  expression  back  into  equation  (132)  for  the  second  derivative  term  and  sim¬ 
plifying  reduces  to 

(kidmil))  dl.  (134) 

Jin 

where  (kidmil))  =  Ji  (kidm(l))  -  jYi  (kidm(l))  is  the  first  order  Hankel  function  of 
the  second  kind. 

Equation  (134)  can  now  be  easily  implemented  in  a  computer  code.  However,  care 
must  be  taken  when  m  =  n,  that  is  when  the  integration  is  performed  over  a  segment 
containing  the  current  match  point.  When  m  =  n,  dm(l)  =  0  at  each  matchpoint 

expression  in  equation  (134)  becomes  a  singiilar  integral.  There 
are  many  techniques  for  handling  the  numerical  integration  of  such  singular  integrals.  One 
of  the  most  efiicient  techniques  is  to  remove  the  singularity  from  the  integrand,  integrate 
the  singularity  analytically  if  possible,  then  add  its  contribution  back  to  the  original  inte¬ 
gral.  In  order  to  apply  this  technique  to  the  current  problem,  the  singular  natmre  of  the 
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integrand  must  first  be  imderstood.  In  the  limit  as  the  argument  of  the  Neumann  function 
goes  to  zero  [6], 

«  —  as  a;  -4-  O'*'.  (135) 

TTX 

Given  this  asymptotic  approximation,  consider  removing  the  singularity  from  the  integrand 
of  equation  (134)  in  the  following  manner. 


<^mn  —  / 

Jin 


^"+*  ki 


Jl{kldmil))-jYUkidm{l))- 


wkidmil)] 


dl 


+  lim 

<5^0+ 


rln+l 

Jin 


(136) 


27r  [dm{l)  +  5] 


Tdl. 


Now  the  first  integral  in  equation  (136)  can  be  integrated  numerically,  and  the  second 
integral  can  be  computed  analytically  to  determine  its  contribution  to  the  entire  integral. 
This  is  done  by  performing  the  integration  and  then  taking  the  limit  as  S  goes  to  0  as 
shown  below. 


lim 


27r  [dm{l)  +  Sf  27r  In+i  In 

-2 _ 2  1  ^  -2 

27r  In+l  ^71  ^n+1  J  ^  (^n+l  ^n) 


(137) 


So  in  order  to  numerically  calculate  the  diagonal  entries  of  the  amn  submatrix  (where 
m  =  n),  the  following  formula  can  be  used 


^mn 


rln+i 

Jin 


ki 


^jdmiP) 


Jl(Mm(0)-J>l(Mm(0)- 


2:/ 


dl- 


1^kidm{l)\  TT  (Zn.f  1  In) 


(138) 


The  Pmn  elements  of  Ai  are  found  using  the  formula 


dydz 


dl.  (139) 


Notice  that  in  order  to  calculate  these  integrals  it  will  be  necessary  to  find  an  expression 
for  the  components  of  the  outward  pointing  normal  vector  along  dD  as  well  as  another 
partial  derivative  of  the  two-dimensional  Green’s  function.  First,  to  obtain  the  normal 
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vector,  define  a  unit  vector  that  is  tangent  to  the  perimeter  of  the  cavity  pointing  in  the 
clockwise  direction.  That  is, 


In 


ivn+l  -yn)y+  (^n+1  -  Zn)  Z 
^n+1  ~  ^n 


(140) 


Then  the  rmit  normal  vector  pointing  out  from  the  surface  and  aperture  of  the  cavity  can 
be  found  by  taking  the  cross  product  of  x  with  1„.  This  results  in 


A  ?  A  ,  Aj/n-fl  Vn 

n  =  X  X  1„  =  -y- - —  +  z- - — — . 

^n+1  m-fl 


(141) 


5j/2 


ihdrail))  +  ihdmil)) 


Furthermore,  it  can  be  shown  that 


{kidmjl))  _ 

dydz 


k\ 


(y(0-ym+i)  {zil)-z^^) 


[dm{l)r 


(kidmil)) 


(142) 


This  gives  the  components  of  the  unit  normal  vector  required  in  equation  (139).  Specifi¬ 
cally,  they  are 

_  Vn+l  ~  Vn 

In+l-ln 

Given  that 

it  was  found  in  the  previous  section  that 

(kidmil)) 


(143) 


(144) 


(145) 
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Substituting  these  expressions  back  into  equation  (139)  and  simplifying  the  results  gives  a 
more  useful  form  of  the  Pmn  equation.  That  is, 


kj  (2) 


r  (kldmil)) 


+  (kidmil)) 


[MDf 

(9(0 -!(„+§)  ^ 


-Uz  +  n. 


[y{l)-ym+i), 


[dm{l)f 


{y{l)-ym+iy 


)  (kMi))  {y{l)  -  ym^)  (^(0  -  dl. 


[dm(0] 


(146) 


Equation  (146)  can  now  be  easily  implemented  in  a  computer  code  to  calculate  the  Prnn 
elements  of  the  impedance  matrix.  Also,  since  m  =  1, 2, ... ,  Ni  and  n  =  iVi  +  1,  iVi  + 
2,...,Ni  +  N2,  pH)  €  S  and  p'{l)  ^  In  other  words,  the  match  points  for  the  ^mn 
entries  all  lie  on  the  cavity  aperture  cr  and  the  integrations  are  all  performed  over  segments 
on  the  cavity  surface  S.  So  the  metric  dm(0  Is  never  0  and  the  integrand  in  equation  (146) 
is  never  singular. 


4-4-^  Computing  A2.  The  first  Ni  by  Ni  submatrix  of  the  A2  matrix  is  composed 
of  all  zeros  except  for  the  diagonal  elements.  They  are  all  constant  and  are  given  by  the 
expression 

H=^{kiZi  +  kQZo),  (147) 

which  can  be  easily  calculated  in  the  computer  code. 

Recall  that  the  rjmn  terms  are  given  by  the  formula 


Vnin  — 


dn 


dl. 


(148) 


Rewriting  this  equation  with  the  equivalent  form  of  the  two-dimensional  Green’s  function 
gives 

_  ^  {kidmji)) 

Jin 


Vmn 


dn 


-dl. 


(149) 
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which  requires  the  computation  of  the  normal  derivative  of  the  Hankel  function.  The 
normal  derivative  of  a  continuous  function  /  with  continuous  first  and  second  partial 
derivatives  in  a  region  defined  by  a  simply  connected  path  C  is  given  by  the  formula 

^-fi-V/,  (150) 

where  n  is  the  outward  pointing  unit  normal  vector  to  C  [20].  In  order  to  compute  the 
normal  derivative  in  equation  (149),  define  the  unit  normal  vector  in  the  same  manner  as 
the  previous  section.  That  is, 


I  ^  yn-\-l  Vn 
In+l-ln' 


(151) 


The  gradient  of  the  Hankel  function  can  be  shown  to  be 

ihdmil))  =  {hdm{l))y 

^  (152) 

Taking  the  dot  product  of  n  with  equation  (152),  plugging  back  into  equation  (149)  and 
simplifying  results  in 


W  =  ^  ^  (fcl<^m(0)  [-%  (y(0  -  ym+i)  -  (^(0  -  ^m+l)]  (153) 

which  is  now  in  a  form  that  can  be  computed  numerically.  Notice  that  in  this  case,  p{l)  €  S 
and  p'{l)  6  a  since  all  of  the  matchpoints  are  contained  in  segments  edong  the  aperture 
while  the  integrations  are  performed  over  segments  along  the  smrface  S.  So  the  expression 
dmii)  will  never  be  zero  and  there  are  no  singularities  to  consider. 

The  7^71  elements  of  the  A2  matrix  eire  given  by  the  formula 


jT'"*'  [klG^  -  fcSGc  (i,i„+i)]  <«. 


Tmn  — 


juei 
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which  when  rewritten  with  the  appropriate  Hankel  functions  substituted  in  for  the  two- 
dimensional  Green’s  functions  appears  as 

Iran  =  ^  f {hdm{l))  -  (Mm(/))1  dl.  (155) 

In  the  case  where  m^n,  equation  (155)  can  be  readily  implemented  in  a  computer 
code  to  produce  the  '^mn  elements  of  A^.  However,  since  p{l)  6  a  and  p'{l)  e  cr,  dm{l) 
will  be  zero  at  the  matchpoints.  This  will  make  the  Neumann  functions  in  the  integrand 
of  equation  (155)  singular  when  m  =  n.  In  order  to  efficiently  calculate  the  numerical 
integration  in  these  instances,  again  consider  removing  the  singularity,  integrating  it  anar 
lytically,  and  adding  its  contribution  back  to  the  total  integral.  Notice,  however,  that  in 
the  degenerate  case  {ko  =  h,  i.e.  an  unfilled  cavity)  the  integrand  in  equation  (155)  will 
be  zero  for  every  value  of  m  and  n.  In  this  case,  the  7mn  submatrix  may  simply  be  set  to 
all  zeros  with  no  treatment  of  the  singularity  necessary. 

The  asymptotic  small  axgmnent  approximation  of  the  zeroth  order  Neumann  function 
is  given  by 

yo(a:)  «  —  Inf^)  as  a; ->  O'*",  (156) 

where  7  ~  1.7810724  is  Euler’s  coefficient  [20].  Since  the  two  Hankel  functions  appearing 
in  equation  (155)  difier  by  only  the  wave  number  fcf  in  their  argument,  a  general  treatment 
of  how  to  remove  the  singularity  where  the  wave  number  is  simply  k  will  be  given  for 
brevity.  Consider  the  equation 

7mn  =  (kdmil))  dl.  (157) 

Jin 

The  singularity  in  the  Neumaim  function  can  be  removed  by  writing 

TT  2 

Jin  ^  2 


dl 

(158) 
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The  first  integral  in  equation  (158)  can  now  be  numerically  calculated  in  the  computer 
code  without  any  problems,  while  the  second  integral  contains  the  removed  singularity.  In 
order  to  integrate  the  singularity  analytically,  consider  splitting  up  the  natural  logarithm 
in  the  following  manner. 

/sing  =  --  In  ^  +  In  {dm{l))  dl  (159) 

^  Jin  2 

By  breaking  the  integral  up,  the  first  term  in  the  integrand  is  simply  a  constant  and  the 
second  term  can  be  integrated  by  carefully  considering  its  nature.  Begin  by  integrating 
the  constant  piece. 

/sing  =  In  (Wl  -  ^n)  -  I  In  (dmil))  dl.  (160) 

The  natural  logarithm  term  in  the  remaining  integral  is  symmetric  about  the  matchpoint, 
so  the  integral  can  be  rewritten  as 

/sine  =  -|  >”  (  y)  ('»+!  -  -  7  1"  j)  (161) 

or  equivalently  as 

9n  /  4^7  r(^n+l— in)/2 

/sing  =  -f  ( Y  j  i^n+l  I  (0  dl.  (162) 

Integrating  the  natural  logarithm  and  simplifying  gives  the  contribution  of  the  removed 
singularity  as 

j  _  2j(ln-f-l  In) 

-'sing  —  - 

If  —  In  is  denoted  as  A„,  the  treated  form  of  equation  (155)  can  be  written  as 


.  (163) 


7mn  —  ~  ^  {l^ldmil))  +  ^  In 


K 

4:uei 


TT 


1-ln 


fco 


4tujei 


(fcodm(O)  +  ^In 


'ykod-mil)  ^  2jAyi 


1-ln 


)]} 


(164) 
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4.4-3  Computing  >13.  Finally,  the  elements  of  the  matrix  must  be  computed. 
The  iprnn  elements  are  given  by  the  equation 


Jin 


dz 


-dl. 


(165) 


This  is  rewritten  with  the  appropriate  two-dimensional  Green’s  function  as 

,  1  dHf\kidm{l))  „ 

= 

The  partial  derivative  of  the  Hankel  function  with  respect  to  ^  is 

(fcldm(O)  • 


(166) 


agf  (fcld^(O)  _  ^m+i) 


dz 


^m(0 


(167) 


So  equation  (165)  can  be  expressed  as 

fcl  /**”+' 


V’n 


= r 


dm(J) 


Hr  {kidm{l))dl, 


(168) 


which  is  now  in  a  form  that  can  be  readily  implemented  in  a  computer  code.  Notice  also 
that  since  p{l)  €  cr  and  p'{l)  6  S  there  is  no  singularity  to  consider  when  performing  the 
numerical  integration. 

The  ^mn  elements  of  the  matrix  were  foimd  to  be  given  by  the  equation 

rln+i  dGoilJrn+O 


Cmn  ~  / 
Jin 

In  this  equation,  Gc  is  defined  to  be 


dn 


-dl. 


(169) 


Gd  (fci ;  p,  p')  =  G  (fci ;  p,  p')  -  G  (fci ;  p  p'i) 


where 


-*!  t 

pi  =  yy-z'^ 


(170) 


(171) 


41 


is  the  reflection  of  the  position  vector  p'{l)  across  the  aperture  of  the  cavity.  For  calcula¬ 
tions  involving  G{ki\p,p[).,a.  new  metric  for  \p{l)  -  p'{l)\  must  be  employed.  Define 


(172) 


to  be  such  a  metric.  Then  taking  the  gradient  of  G  (fei;  p,  p')  and  of  G  {ki;p,  p'^)  as  before, 
and  using  the  same  unit  normal  vector  as  defined  in  equation  (151),  the  normal  derivative 
of  Gd  (fci;  p,  p\)  can  be  found.  Doing  so  allows  equation  (169)  to  be  rewritten  as 


Cmn 


{kidmil)) 


n. 


(j/CO-ym+j)  (»(0  +  ^mti) 


+  7^- 


d\n{l) 


{ki4n{l))  dl 


(173) 


Even  though  p{l)  G  S  and  p'{l)  6  S,  notice  that  the  newly  defined  metric  in  equation 
(172)  will  never  be  zero.  So  there  will  be  no  trouble  with  singularities  when  numerically 
integrating  the  Hankel  function  that  involves  However,  the  Hankel  function  that 

contains  dm{l)  may  pose  some  problems  when  the  numerical  integration  is  performed  over 
a  segment  contedning  the  current  matchpoint.  So  if  m  =  n,  the  singularity  should  be  ex¬ 
tracted  and  integrated  separately.  The  same  procedme  presented  in  equations  (135)-(138) 
can  be  followed  to  remove  the  singiilarity  by  subtracting  out  the  asymptotic  small  argument 
approximation  of  the  first  order  Neumann  function  from  the  integrand  and  integrating  it 
separately. 

It  should  be  pointed  out  here  that  the  singularities  that  arise  in  filling  the  diagonal 
elements  of  the  amn  and  ^rnn  submatrices  axe  quite  strong.  There  is  another  much  more 
easily  implemented  method  of  dealing  with  such  strong  singularities  that  is  commonly  used 
in  computational  electromagnetics.  Essentially  this  involves  redefining  the  coordinates  of 
the  matchpoint  when  Ccilculating  the  self-impedance  terms  (where  m  =  n).  By  “bmnping” 
the  matchpoint  a  distance  of  say  one  hundredth  of  a  wavelength  away  from  the  cavity 
jdong  the  outward  pointing  normal  vector,  the  metric  dm(0  will  never  actuedly  be  zero. 
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This  will  allow  the  numerical  integration  to  be  performed,  and  any  error  introduced  in  the 
solution  when  the  system  is  solved  should  be  negligible. 

Another  similar  method  of  dealing  with  strong  singular  integrals  in  computational 
electromagnetics  is  the  method  of  equivalent  separation.  Using  this  method,  a  small  in¬ 
cremental  distance  is  added  to  dmil)  on  every  segment  [10].  This  is  another  way  of  not 
allowing  the  argument  of  the  Hankel  functions  to  ever  be  zero,  thereby  obtaining  finite  val¬ 
ues  when  the  integrations  are  performed.  However,  this  technique  introduces  a  phase  error 
into  the  integrands,  which  can  have  a  significant  effect  on  the  accuracy  of  the  solution. 


4-5  Error  Estimate 

In  this  section  an  estimate  of  the  order  of  convergence  of  the  MoM  algorithm  is 
developed.  Consider  using  a  mesh  of  equally  sized  segments  to  discretize  the  perimeter  of 
the  cavity  dD.  Let  h  be  the  length  of  each  segment.  In  the  L^-norm,  the  total  error  in  the 
approximation  of  u{l)  is  given  by 


= 


\u{l)  -  dl. 


(174) 


where  is  the  approximate  average  value  (a  constant)  of  u{l)  on  each  segment  defined  by 
[In,  In+i]  [5]-  Notice  that  the  approximation  as  well  as  the  interval  [/«)  ^n+i]  depend  on 
the  mesh  size  h.  Since  a  point  matching  scheme  was  used  at  the  midpoint  of  each  segment, 
u{cn)  —  Ufi^  for  every  midpoint  c„  =  (/n+i  “  in)/2.  Furthermore,  by  the  Mean- Value 
theorem,  there  exists  a  point  d„  in  each  interval  [Z„,  Z„+i]  such  that 


||eft|U2<  f^^'W'idn)-{l-Cn)\ 


dl. 


(175) 


Taking  the  maximum  of  |/  —  Cn|  on  each  segment  gives 


Wn+1  I  h 


dl. 


(176) 
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Now  taking  the  maximiun  of  u'{l)  on  each  segment  and  doing  the  integration  reduces  the 
inequality  in  (176)  to  _ 

||eft||i2  <  J  XI  “'(0  •  y-  (177) 

V  nedD 

Since  equally  sized  segments  were  used  to  discretize  the  cavity,  the  summation  under  the 
square  root  can  be  replaced  giving 

||eft||£2  <  •  max  u'{l)  •  y,  (178) 

where  F  is  the  total  arc  length  along  the  perimeter  dD.  Finally,  simplifying  gives  the  error 
estimate  in  the  L^-norm  as 

<  I  Vr  •  max  u'{l)  =  C-h^  0{h).  (179) 

According  to  (179),  if  a  total  of  N  segments  are  used  along  dD  to  approximate  u{l),  and 
then  the  mesh  is  refined  by  cutting  the  length  of  each  segment  in  half  for  a  total  of  2N 
segments,  the  error  in  the  approximation  of  u(/)  should  be  cut  approximately  in  half.  This 
estimate  of  the  order  of  convergence  will  be  examined  in  the  next  chapter  by  considering 
the  relative  error  in  the  solutions  of  several  test  cavities. 
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V.  Numerical  Results 


In  this  chapter,  the  integral  equation  method  is  applied  to  several  test  cavities  to  calculate 
their  monostatic  RCS  signatures.  The  results  found  are  then  compared  with  those  foimd 
by  RAM2D  [7]  and  by  a  hybrid  finite  element/Fourier  transform  code  [25],  and  the  method 
is  shown  to  be  very  accurate.  The  chapter  also  contains  some  error  history  plots  which 
agree  with  the  predicted  convergence  estimate  presented  in  the  previous  chapter. 

5.1  Computing  the  Monostatic  RCS 

Once  the  imknown  fields  have  been  computed  for  a  given  incident  angle,  a  method  for 
calculating  the  backscatter  RCS  is  required.  In  order  to  develop  such  a  method,  consider 
the  following  mathematical  definition  of  the  RCS  of  a  two-dimensional  target. 

i#(p')r 

0-2-D  =  ,  fim  27r  |p'|  J - Lg,  (180) 

where  |p'|  =  p'  is  the  distance  from  the  target  to  the  observer,  and  E®  and  E*"*^  are  the 
scattered  and  incident  electric  fields,  respectively  [6].  Recall  the  radiation  integral  given 
by  Application  5  in  equation  (33)  was 

^E^(p')  =  jT  z  X  E(^  •  F(feo;  p, p')da.  (181) 

This  will  give  the  necessary  numerator  in  equation  (180).  The  denominator  is  already 
known  since  it  is  simply  the  given  incident  electric  field.  In  order  to  compute  E®(p')  in 
terms  of  the  magnetic  fields,  apply  Maxwell’s  equation  (2)  to  obteiin 

E®(/)  =  ^  fix  ( —V  X  H(^)  .  f  (fco;  p,  p')d(r.  (182) 

*0  J(r  VJWeo  / 

Letting  H(p)  =  ux,  perform  the  two  vector  cross  products  and  substitute  the  appropriate 
expression  in  for  F  to  reduce  the  equation  to 
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Cross  multiplying  the  gradient  of  the  Green’s  function  with  the  idemfactor  and  dotting 
the  result  with  x  simplifies  the  scattered  electric  field  to 


da. 


(184) 


Denote 

k=^/{y-y'Y  +  {z-z'Y  (185) 


so  that 

GQ  =  jH'i\kQR).  (186) 

Then  taking  the  appropriate  partial  derivatives  of  the  Hankel  function  allows  equation 
(184)  to  be  rewritten  as 


(187) 


Considering  the  geometry  in  Figure  4,  it  is  apparent  that  R  diS  I  p  I  ^  OO  SU1C0 
y  and  z  are  restricted  to  the  aperture  of  the  cavity.  Furthermore,  since  z  G  a,  z  =  0  and 
R  can  be  rewritten  as 

(188) 
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This  can  be  further  simplified  to 


R=  i/y2  _  2y'y  +  (y')2  +  (^')2 


=  y/Wf+W-Wy 

=  V (p')^  +  2/^  -  2yp'  cos(0) 


(189) 


=  p'\  1  + 


a)- 


2y  cos(0) 

~7“ 


=  p'-ycos(0)  +  (!l(jr) 


In  the  limit  as  ^  >  oo,  the  large  argmnent  approximation  of  the  Hankel  function  is 


lim  ~  \ 

R->oo  V 


-jkoR 


(190) 


This  allows  the  scattered  electric  field  to  be  approximated  as 


vS /  (-y^  - (191) 


Recall  the  incident  electric  field  was  given  as 


W°%p')  =  Zo(sin0y  -  cos0z)e^*»(»®’"®+^'=“®''). 


(192) 


Now,  plugging  equations  (191)  and  (192)  into  (180)  and  simplifying  yields 


<72.D  ~  27r 


1  f  ^^koycos{0)  jj 

°  V  dz  ^  ^R3/2 


.  -  ^/]^(y  -  y') 


R3/2 


) 


rfCT 


|Zo(sin0y  —  cos  0z)| 


(193) 


Recognizing  that  lim  J^/y/k  =  1,  and  lim  y/R  =  0,  and  again  applying  the  geometry 

P'-^OO 


p'->CX) 


in  Figure  4,  this  can  be  rewritten  as 

sin(0)y  +  cos(0)z)  ^  |^e^*oycos(0)^^ 


0’2-D  ~  27r 


|(sin0y  —  cos  0z)|'‘ 


(194) 
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Finally,  converting  the  integral  into  a  summation  and  simplifying  the  result  reduces  the 
RCS  equation  to 


1 

o'2-D  ~  7- 
ko 


n=l 


(195) 


where  the  bn  are  the  approximations  to  ^  on  each  segment,  A„  =  In+i  —  In  is  the  length 
of  each  segment,  is  the  matchpoint  of  each  segment,  and  0  is  the  angle  of  incidence. 


5.2  RCS  Plots 

The  integral  equation  method  is  applied  to  several  test  cavities  in  this  section,  and  the 
monostatic  RCS  signatures  axe  computed  for  incident  angles  between  90°  and  180°.  The  ac¬ 
curacy  of  the  method  is  considered  by  comparing  the  plots  with  the  RCS  values  calciilated 
by  RAM2D  [7]  as  well  as  with  the  results  from  a  hybrid  finite  element/Fourier  transform 
code  written  by  Van  [25].  The  RAM2D  code  is  a  commercial  grade  two-dimensional  scat¬ 
tering  problem  solver,  and  Van’s  finite  element  code  has  been  benchmarked  against  other 
existing  known  results. 

5.2.1  Test  Cavity  1.  For  the  first  test  case,  consider  the  cavity  in  Figure  5.  It 
is  a  rectangular  cavity  with  dimensions  1  meter  wide  by  0.25  meters  deep.  The  interior 
of  the  cavity  is  unfilled,  which  is  denoted  by  the  free  space  parameters  (eo,/2o)-  So  in  this 
case,  ei  =  eo  and  ni  =  fio-  Let  an  incident  plane  wave  imder  TE  polarization  and  with 
frequency  300  MHz  impinge  on  the  cavity  to  produce  the  scattered  fields.  The  incident 
wave  will  have  a  wavelength  of  A  =  1  meter,  emd  its  wave  number  or  propagation  constant 
will  be  fc  =  27r.  Notice  that  since  the  cavity  is  unfilled,  the  incident  field’s  propagation 
constant  is  the  same  in  the  upper  half  plane  as  it  is  within  the  interior  of  the  cavity. 

The  MoM  algorithm  was  employed  using  10  pulse  basis  functions  per  wavelength. 
This  results  in  having  10  segments  on  the  apertmre  {Ni  =  10)  and  16  segments  on  the 
surface  of  the  cavity  (IV2  =  16),  making  =  /  a  36  by  36  linear  system.  The  unknown 
expansion  coefficients  were  approximated  by  the  method  to  produce  the  total  magnetic 
field  on  the  perimeter  of  the  cavity  and  the  normal  derivative  of  the  total  magnetic  field  on 
the  aperture  of  the  cavity.  The  magnitudes  of  these  field  quantities  (in  Amperes/meter) 
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are  shown  in  Figures  6,  7,  and  8  for  an  incident  angle  of  0  =  90°.  (It  should  be  pointed  out 
here  that  more  than  10  pulse  basis  functions  per  wavelength  were  used  to  produce  the  plots 
of  the  fields  in  the  following  3  figmes.  Actually,  25  pulses  were  used  so  that  the  graphs 
would  be  fairly  smooth.  However,  only  10  pulses  per  wavelength  were  used  to  compute  the 
RCS  values.) 

In  Figures  6,  7,  and  8  notice  that  the  fields  all  have  a  symmetric  shape  with  respect 
to  the  center  of  the  cavity.  This  is  to  be  expected  when  a  symmetric  cavity  is  illuminated 
by  a  plane  wave  at  normal  incidence.  Finally,  in  Figure  9  the  graph  of  the  monostatic  RCS 
values  computed  by  the  integral  equation  method  is  plotted  for  several  angles  between 
90°  and  180°.  The  integral  equation  method  output  is  compared  with  the  RCS  values 
computed  by  RAM2D  in  the  figure,  and  excellent  agreement  is  observed.  So  the  MoM 
algorithm  was  able  to  produce  very  accurate  results  while  only  having  to  solve  a  36  by  36 
linear  system  for  each  incident  angle. 
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OO 


Magnitude  of  total  magnetic  field  on  aperture 


Figiire  6  Magnitude  of  the  total  magnetic  field  on  the  aperture,  normal  incidence 
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Cavity  surface,  S 

Figure  7  Magnitude  of  the  total  magnetic  field  on  the  cavity  surface,  normal  incidence 
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Magnitude  of  normal  derivative  of  magnetic  field  on  aperture 


Figure  8  Magnitude  of  the  normal  derivative  of  the  magnetic  field  on  the  aperture, 
normal  incidence 
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Monostatic  RCS  (unfilled  rectangular  cavity) 


Figure  9  Monostatic  RCS  signature  of  an  unfilled  rectangular  cavity.  Integral  equation 
method  (o)  and  RAM2D  ( — ) 
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5.2.S  Test  Cavity  2.  For  the  second  test  case,  another  unfilled  rectangular  cavity 
is  illuminated  by  a  300  MHz  plane  wave  tmder  TE  polarization.  Consequently,  the  wave¬ 
length  of  the  incident  field  is  A  =  1  meter  and  the  propagation  constant  is  again  k  =  2n 
both  in  the  upper  half  plane  and  the  interior  of  the  cavity.  However,  this  time  the  cavity 
is  much  deeper  and  a  little  bit  wider.  As  shown  in  Figiure  10,  the  cavity  is  1.2  meters  wide 
and  0.8  meters  deep. 


The  MoM  algorithm  was  again  employed  using  10  pulse  basis  functions  per  wave¬ 
length.  Due  to  the  geometry  of  the  cavity,  this  resulted  in  a  total  of  12  segments  on  the 
cavity  aperture  (iVi  =  12)  £ind  28  segments  on  the  cavity  surface  {N2  =  28).  So  the  total 
nximber  of  unknown  expansion  coefiicients  is  52,  making  Au  =  /  a  52  by  52  linear  system. 
Again,  using  equation  (195)  the  monostatic  RCS  values  were  computed  by  the  MoM  al¬ 
gorithm  for  incident  angles  between  90°  and  180°,  and  the  results  were  compared  to  the 
RCS  values  computed  by  RAM2D.  In  Figure  11,  the  two  sets  of  data  are  compared,  and 
excellent  agreement  between  the  two  methods  is  again  observed. 
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Monostatic  RCS  (deep  unfilled  rectangular  cavity) 


Figiire  11  Monostatic  RCS  signature  of  a  deep  unfilled  rectangular  cavity.  Integral  equa¬ 
tion  method  (o)  and  RAM2D  { — ) 


55 


5.2.3  Test  Cavity  3.  For  the  third  test  case,  a  rectangular  cavity  filled  with 
a  homogeneous,  lossy  material  is  considered.  The  incident  field  for  the  problem  is  still 
a  TE  polarized  plane  wave  with  frequency  300  MHz  and  wavelength  A  =  1  meter.  The 
propagation  constant  in  the  upper  half  plane  is  still  fco  =  27r,  however  the  fields  will  have  a 
different  wave  number  within  the  cavity  due  to  the  presence  of  the  lossy  matericd.  In  this 
case,  the  cavity  is  filled  with  a  lossy  dielectric  material  characterized  by  er  =  4  —  y,  so  the 
propagation  constant  within  the  material  filled  cavity  is  fci  =  27rv'4  —  j.  The  dimensions 
of  the  cavity  are  the  same  as  the  rectangular  cavity  from  the  first  test  case,  1  meter  wide 
and  0.25  meters  deep.  The  geometry  of  the  problem  is  shown  in  Figmre  12. 


The  monostatic  RCS  values  were  computed  and  compared  with  the  results  found  by 
a  hybrid  finite  element/Fourier  transform  code  written  by  Van  [25].  The  plots  appear  in 
Figure  13  and  indicate  very  good  agreement  between  the  two  methods.  In  comparing  the 
two  RCS  signatures  in  Figures  9  and  13,  the  effect  of  the  lossy  dielectric  material  appears 
to  be  fairly  significant.  Both  graphs  have  a  null  appearing  around  9  =  120°  and  have  very 
similar  overall  shapes,  but  the  RCS  values  for  the  filled  case  are  much  lower  than  for  the 
unfilled  case.  This  is  because  the  dielectric  material  absorbs  some  of  the  energy  of  the 
incident  field  which  results  in  a  decreased  scattered  field. 
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Monostatic  RCS  (filled  rectangular  cavity) 


Figure  13  Monostatic  RCS  signature  of  a  rectangular  cavity  filled  with  a  lossy  dielectric 
material.  Integral  equation  method  (o)  and  finite  element  method  ( — ) 
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5.24  Test  Cavity  4.  The  fourth  test  cavity  is  a  deep  rectangular  cavity  filled 
with  a  homogeneous,  lossy  dielectric  material.  Again,  the  incident  field  is  a  TE  polarized 
plane  wave  with  frequency  300  MHz.  The  relative  electric  permittivity  of  the  medium  is 
Cr  =  2  —  j  which  makes  the  wave  number  ki  =  2ir\/2  —  j  inside  the  cavity  interior.  In  the 
upper  half  plane,  the  propagation  constant  is  stUl  feo  =  27r,  and  the  incident  field  has  a 
wavelength  of  A  =  1  meter.  The  geometry  of  the  fourth  test  case  is  shown  in  Figure  14. 
The  cavity  is  1.2  meters  wide  and  0.8  meters  deep. 


The  MoM  algorithm  was  employed  using  10  pulse  basis  functions  per  wavelength. 
This  resulted  in  solving  a  52  by  52  linear  system  for  each  incident  angle.  The  RCS  val¬ 
ues  given  by  the  integral  equation  method  as  well  as  a  comparison  with  Van’s  finite  ele¬ 
ment/Fourier  transform  method  are  shown  in  Figure  15.  Again,  for  a  very  small  compu¬ 
tational  effort,  strong  agreement  is  observed  between  the  two  methods. 
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Monostatic  RCS  (deep  filled  rectangular  cavity) 


Figure  15  Monostatic  RCS  signature  of  a  deep  rectangrilar  cavity  filled  with  a  dielectric 
material.  Integral  equation  method  (o)  and  finite  element  method  ( — ) 
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5.2.5  Test  Cavity  5.  The  final  test  cavity  is  a  triangular  shaped  cavity  filled  with 
a  homogeneous,  lossy  dielectric  material.  The  dimensions  of  the  triangular  cavity  are  1 
meter  wide  across  the  aperture  and  0.5  meters  deep.  The  interior  of  the  cavity  is  filled 
with  a  dielectric  material  having  relative  electric  permittivity  er  =  4  —  j.  This  makes  the 
wave  number  in  the  interior  of  the  cavity  ki  =  2ir —  j.  The  incident  field  is  a  300  MHz 
plane  wave  under  TE  polarization.  Agciin,  the  wavelength  of  the  incident  field  is  A  =  1 
meter,  and  its  free  space  propagation  constant  is  ko  =  2iv.  The  geometry  of  the  fifth  test 
case  appears  in  Figure  16. 


Using  10  pulse  basis  functions  per  wavelength,  the  perimeter  of  the  cavity  is  dis¬ 
cretized  into  10  segments  along  the  aperture  and  16  segments  along  the  surface.  This 
results  in  the  MoM  algorithm  solving  a  36  by  36  linear  system  for  each  incident  angle. 
Comparing  the  RCS  values  of  the  integral  equation  method  with  those  found  by  the  finite 
element /Fomier  transform  method  in  Figure  17  again  shows  very  good  agreement. 
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RCS  (dbm) 


Monostafic  RCS  (filled  triangular  cavity) 


Figure  17  Monostatic  RCS  signature  of  a  triangular  cavity  filled  with  a  dielectric  mate¬ 
rial.  Integral  equation  method  (o)  and  finite  element  method  ( — ) 
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5.3  Error  History  Plots 

In  this  section,  the  rate  of  convergence  of  the  integral  equation  method  is  examined. 
In  the  following  experiments,  each  cavity  is  initially  meshed  with  N  equally  sized  segments 
of  length  h,  and  the  MoM  algorithm  is  used  to  compute  the  total  magnetic  field  along  the 
cavity  boundary  dD.  Then,  the  mesh  is  refined  by  chopping  each  segment  in  half  so  that 
there  are  a  total  of  2N  segments  along  dD,  and  the  total  magnetic  field  is  recomputed.  At 
each  step  the  relative  error  in  the  L^-norm  is  computed  using  the  formula 


error^  = 


\\Uhn  -  Uhn-liy 

W'^hnWl^ 


(196) 


where  hn  =  hn_i/2  indicating  that  the  mesh  is  chopped  in  half  at  each  step,  is  the 
approximate  solution  at  the  current  size  step,  and  is  the  approximate  solution  at  the 
previous  size  step. 

Error  history  -  Normal  liciderKO 


Figure  18  Relative  error  history  of  an  empty  rectangular  cavity  for  normal  incidence 


Test  Case  1.  For  the  first  experiment,  consider  the  scattering  fi-om  an  unfilled  cavity  1 
meter  wide  by  0.25  meters  deep  of  a  300  MHz  plane  wave  imder  TE  polarization.  Let  the 
plane  wave  impinge  on  the  cavity  imder  normal  incidence  {9  =  90°).  Using  equation  (196), 
the  relative  error  history  versus  1/h  is  plotted  in  Figure  18,  and  a  log2-log2  plot  of  the 
relative  error  versus  1/h  is  given  in  Figure  19. 
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Error  Plot:  Normal  tnctdanc©  (Log-Log) 


Figure  19  Log-Log  plot  of  relative  error  history  for  normal  incidence 

Notice  in  Figure  18  that  at  each  size  step  where  h  is  halved,  the  relative  error  seems 
to  be  cut  approximately  in  half.  So  for  h  sufficiently  small  the  L"^  error  rate  is  approaching 
1  for  the  solution  of  an  unfilled  rectangular  cavity  under  normal  incidence.  Furthermore, 
the  slope  of  the  line  in  Figure  19  is  approximately  -1.07  which  is  indicative  of  an  0{h) 
method.  This  agrees  with  the  estimate  given  in  the  previous  chapter. 


Error  History:  10  degreos 


Figure  20  Relative  error  history  of  an  empty  rectangular  cavity,  0  =  10° 

Test  Case  2.  In  this  experiment,  consider  the  same  cavity  as  in  Test  Case  1  with  a  300 
MHz  plane  wave  under  TE  polarization  impinging  on  the  cavity  at  an  incident  angle  of 
0  =  10°.  According  to  Figure  20,  the  error  rate  also  appears  to  be  approaching  1  for 
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this  case  as  the  mesh  is  refined  at  each  size  step.  And  the  slope  of  the  line  in  Figure  21  is 
approximately  -1.1  which  again  agrees  with  the  rate  of  convergence  estimate. 


Error  Plot:  10  degrees  (Log-Log) 


Figure  21  Log-Log  plot  of  relative  error  history,  9  =  10° 


Error  History  for  Filled  Triangular  Cavity 


Figure  22  Relative  error  history  of  a  filled  triangular  cavity,  6  =  35° 

Test  Case  3.  In  this  experiment,  a  filled  triangular  cavity  is  considered.  The  cavity  is  in 
the  shape  of  an  isosceles  triangle  1.2  meters  by  0.8  meters  deep.  It  is  filled  with  a  dielectric 
material  having  relative  electric  permittivity  =  2.26.  Again  the  incident  wave  is  a  TE 
polarized  plane  wave  with  frequency  300  MHz,  £ind  it  is  propagating  at  an  incident  angle 
of  0  =  35°.  The  relative  error  plot  in  Figme  22  again  indicates  that  for  h  sufiiciently  small, 
the  rate  of  convergence  is  approaching  1  since  the  relative  error  is  cut  approximately  in 
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half  at  each  step.  Also,  the  log2-log2  plot  of  the  relative  error  versus  l//i  in  Figure  23 
shows  a  slope  of  about  -1.02  which  is  consistent  with  an  0(h)  rate  of  convergence. 


Enof  Ptot:  Filled  Triangutar  Cavity  (Log-Log) 


Figure  23  Log-Log  plot  of  relative  error  for  a  filled  triangular  cavity,  0  =  35° 


65 


VI.  Conclusions 


This  research  has  investigated  the  application  of  a  moment  method  based  integral  equation 
scheme  to  compute  the  scattering  of  a  TE  polarized  plane  wave  from  a  cavity  embedded 
in  a  PEC  ground  plane.  This  is  motivated  primarily  by  the  need  of  the  United  States 
military  to  accurately  model  and  predict  the  RCS  signatures  of  targets.  The  reduction 
and/or  enhancement  of  a  friendly  air  vehicle’s  RCS  signatmre  as  well  as  the  detection  and 
identification  of  an  enemy  air  vehicle  are  crucied  to  the  success  of  an  air  campaign.  Cavities 
on  an  aircraft,  such  as  the  jet  engine  inlet  ducts  and  exhaust  nozzles,  can  dominate  its 
RCS  signatme.  This  coupled  with  the  challenging  computational  nature  of  cavities  make 
the  study  of  the  scattering  from  cavities  vitally  important. 

There  is  a  large  body  of  research  regarding  the  cavity  scattering  problem,  and  there 
are  many  good  methods  for  computing  the  RCS  of  cavities.  Among  these  are  the  gen¬ 
eralized  network  formulation  (CNF),  integral  equation  methods  such  as  the  electric  field 
integral  equation  (EFIE)  method  and  the  magnetic  field  integral  equation  (MFIE)  method, 
and  the  hybrid  finite  element  methods.  One  of  the  primary  advantages  of  the  method  pre¬ 
sented  here  is  that  it  does  not  suffer  from  the  problem  of  spurious  resonances  at  certain 
frequencies  like  some  of  the  other  methods.  Also,  it  is  very  computationally  efficient  since 
only  the  boundary  of  the  cavity  must  be  meshed  to  solve  the  problem.  Once  the  fields 
along  the  boimdary  of  the  cavity  are  found,  the  fields  everywhere  above  the  cavity  can 
be  computed  for  applications  such  as  capturing  the  RCS  signature  of  the  cavity.  So  very 
accurate  results  can  be  produced  with  a  minimal  amount  of  computations  required. 

In  the  previous  chapters,  the  fundamental  theorems  developed  in  [4]  and  [26]  were 
converted  to  vector  integral  equations  that  govern  the  electromagnetic  scattering  from  a 
cavity  embedded  in  a  PEC  grormd  plane.  A  coupled  set  of  scalar  integral  equations  were 
then  derived  for  the  problem  of  a  TE  polarized  plane  wave  impinging  on  a  two-dimensional 
cavity,  where  the  cavity  could  either  be  empty  of  filled  with  a  homogeneous  material.  These 
scalar  integral  equations  were  then  solved  via  a  Method  of  Moments  algorithm,  and  the 
RCS  signatures  of  several  test  cavities  were  computed.  The  results  produced  compcired 
very  well  with  the  results  of  other  methods  including  RAM2D  and  a  finite  element  method 
code,  and  in  doing  so,  the  method  was  very  computationally  eflftcient.  An  estimate  of  the 
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order  of  convergence  of  the  method  was  also  developed,  and  the  nmnerical  results  indicated 
that  it  was  an  accurate  estimate. 

This  research  represents  a  follow-on  effort  to  work  previously  accomplished  by  Asves- 
tas,  Kleinman,  Wood,  and  Wood  [4],  [26],  [27],  but  the  problem  certainly  has  not  been 
exhausted.  There  are  several  areas  where  future  research  could  be  accomplished.  First, 
instead  of  simply  having  an  empty  or  filled  cavity,  a  thin  dielectric  coating  on  the  surface 
of  the  cavity  could  be  modeled  or  a  cavity  filled  with  layers  of  different  materials.  Also, 
another  area  of  research  would  be  in  comparing  the  rate  of  convergence  of  the  method  with 
the  rates  of  other  methods.  Furthermore,  while  the  method  has  been  shown  to  produce 
accurate  results,  it  would  be  useful  to  prove  that  the  method  always  obtains  the  unique 
solution  to  the  scattering  problem.  Finally,  the  problem  considered  in  this  research  could 
be  extended  to  the  three-dimensional  case,  since  the  integral  equations  are  valid  in  both 
two  and  three  dimensions. 
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